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ABSTRACT  * 


This  report  presents  general  procedures  for  solving  field  problems 
of  engineering  interest  using  digital  computer  techniques.  The  basic 
concept  is  to  represent  a  boundary-value  problem  by  a  superposition 
integral,  approximate  the  integral  equation  by  a  matrix  equation,  and 
invert  the  matrix  for  a  solution.  The  theory  is  dc  scribed  in  terms  of 
the  method  of  moments,  which  is  equivalent  to  the  variational  method. 

For  electromagnetic  antenna  and  scattering  problems,  the  method 
gives  a  matrix  whose-  elements  can  be  interpreted  as  generalized  impedances. 
These  impedances  are  closely  related  to  those  used  in  the  theory  of 
loaded  antennas  and  seatterers,  and  hence  such  loaded  structures  can  also 
be  treated.  A  solution  for  wire  antennas  and  seatterers  of  arbitrary 
shape  is  formulated  in  detail,  and  calculations  for  linear  wire  antennas 
and  seatterers,  both  loaded  and  unloaded,  have  been  made. 

Additional  problems  treated  by  these  procedures  are  two-dimensional 
scattering  by  conducting  cylinders  and  by  dielectric  cylinders,  end 
three-dimensional  scattering  by  bodies  of  revolution.  These  problems 
are  used  to  show  the  effect  of  various  approximations  in  the  solution, 
in  an  attempt  to  draw  some  general  conclusions  as  to  the  best  approxi¬ 
mations.  Special  procedures  for  inverting  matrices  have  also  been  considered, 
to  take  into  account  any  symmetry  properties  present  in  the  matrices. 

A  considerable  saving  in  computation  time  can  often  be  made  by  properly 
utilizing  these  symmetries. 
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The  objective  of  this  effort  was  to  develoe  trace! cal  tec hni cues 


'or  solving  a  wise  class  of  electromagnetic  seaiterli 


s-s» 
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A  unified  ana  general  approach  is  takes  to  electromagnetic  field 
problems  is  terms  of  numerical  analysis.  The  approach  results  Is 
matrix  equations  representing  the  solution  of  the  integro-differsstial 
equations  of  Hie  problem  being  analyzed.  The  accuracy  of  the 
solution  obtained  is  very  good  for  most  engineering  purposes.  Is 
some  cases ,  specifically  that  case  of  the  li near  wire  antenna,  it 
appears  that  the  solution  is  mure  accurate  than  any  previously 
obtained.  The  unifies  viewpoint  developed  is  this  effort  con¬ 
tributes  additional  insight  into  the  solution  of  elec tronagneti c 
scattering  process,  ss  establishes  a  general  framework  with  which 


tc  conduct  the  analysis.  Techniques  and  ecus  applications  of  the 
method  of  approach  are  described  in  Volnse  I  of  the  final  resort. 

The  results  of  this  effort  are  directly  applicable  to  the 
problem  of  the  linear  wire  antenna,  in  that  computer  programs  ass 
complete  numerical  results  are  presented  is  Volume  II  of  the  final 
resort. 
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MATRIX  METHODS  FOR  SOLVING  FIELD  PROBLEMS 

I .  BASIC  CONCEPTS 

A.  Introduction.  The  object  of  this  project^ is  to  develop  practical 

"  *  ^ 

techniques  for  solving  electromagnetic  scattering  problems  using  digital 
computers.  The  general  approach  is  to  approximate  the  integro -differential 
equations  describing  the  boundary  value  problem  by  matrix  equations  and 
to  invert  the  matrix  equations. 

The  theory  is  conveniently  discussed  using  the  concepts  of  linear 
spaces.  The  equations  of  interest  are  of  the  general  type 

L(f)  =  g  (1-1) 

where  L  is  a  known  linear  operator,  g  is  a  known  function,  and  f  is  the 
unknown  to  be  determined.  The  general  procedure  will  be  to  take  an 
equation  involving  functions,  such  as  (l-l),  and  to  approximate  it  by 
a  matrix  equation.  The  solution  to  an  inhomogeneous  matrix  equation  is 
found  by  matrix  inversion.  Most  computer  libraries  have  subroutines  for 
matrix  inversion,  and  hence  a  problem  will  be  considered  solved  once  the 
elements  of  a  well  conditioned  matrix  are  evaluated. 

B.  Formulation  of  Problems.  General  procedures  for  solutions  will 
be  discussed  in  the  language  of  linear  spaces  and  operators.  Hence,  the 
problems  to  be  considered  should  be  formulated  in  this  notation.  Given 
a  deterministic  problem,  we  wish  to  put  it  into  the  form  L(f)  =  g, 
identifying  the  operator  L,  its  domain  (the  functions  f  on  which  it 
operates),  and  its  range  (the  functions  g  resulting  from  the  operation). 
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Furthermore,  we  usually  need  an  Inner  product  <f,g>,  defined  to  satisfy^ 


<f,g>  =  <g,f> 

(1-2) 

<0ff  4  pg,h>  =  Of<f  ;h>  4 

0<g,h> 

(1-3) 

<f*,f>  >  0 

if  f  ^  0 

=  0 

if  f  =  0 

(1-4) 

where  a  and  p  are  scalars,  and  *  denotes  complex  conjugate.  We  some- 

3. 

times  need  the  adjoint  operator  L  and  its  domain,  defined  by 

<Lf,g>  =  <f,Lag>  (1-5) 

Si 

for  all  f  in  the  domain  of  L.  An  operator  is  self  adjoint  if  L  =  L  and 
the  domain  of  I.  is  that  of  L. 

Properties  of  the  solution  depend  upon  properties  of  the  operator. 

An  operator  is  real  if  Lf  is  real  whenever  f  is  real.  An  operator  is 
positive  definite  if 

<f*,Lf>  >  0  (1-6} 

for  all  f  /  0  in  its  domain.  It  is  positive  semidefinite  if  >  is  replaced 
by  >  in  (1-6),  negative  definite  if  >  is  replaced  by  <  in  (.1-6),  etc. 

We  shall  identify  other  properties  of  operators  as  we  need  them. 

If  the  solution  to  L(f)  =  g  exists  ana  is  unique  for  all  g,  then 
the  inverse  operator  L  1  exists  such  that 

f  =  L_1(g)  (1-7) 

If  g  is  known,  then  (l-7)  represents  the  solution  to  the  original  problem. 
However,  (1-7)  is  itself  an  inhomogeneous  equation  for  g  if  f  is  known, 
and  its  solution  is  L(f)  =  g. 

The  usual  definition  of  inner  product  in  Hilbert  space  corresponds  to 
<f*,g>  in  our  notation.  For  this  report  it  is  more  convenient  to  show 
the  conjugate  operation  explicitly  wherever  it  occurs,  and  to  define 
the  adjoint  operator  without  conjugation. 
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In  general,  it  is  desirable  to  think  of  L  and  L-1  as  a  pair  of  operators, 
each  one  of  which  is  the  inverse  of  the  other.  The  choice  in  any  parti¬ 
cular  problem  is  one  of  convenience. 


C.  Example;  Electrostatics .  Consider  a  volume  density  of  charge 
p(x,y,z)  in  unbounded  Euclidean  space  of  constant  permittivity  e,  and 
its  associated  electrostatic  potential  0{x,y,z).  The  differential 
equation  for  the  problem  is  Poisson’s  equation 

-e  ^0  =  p  (1-8) 

subject  to  the  boundary  condition  that  r0-*C1  as  r— ♦<»  for  every  p  of 

finite  extent  and  where  C  ^  a  A 

2  is  a  constant.  In  operator  notation,  ( 1—8)  is 

L0  =  p  (1-9) 

where  „ 

X,  =  -€  v  (1-10) 


subject  to  the  above  boundary  condition.  The  well-known  solution  to 
(1-8)  is 


0(x,y,z) 


i  \ 
/ 


dx*  dy*  az* 


(l-Xl) 


//  2  2  '2 

where  R  =  j(x-x')  +  (y-y*)  +  (z-z*)  is  the  distance  from  &  source 

point  (x*,y*,z*)  to  a  field  point  (x,y,z).  Hence,  the  inverse  operator 


is 


dy*  dz  * 


a-i2) 


It  is  important  to  realize  that  (1-12)  is  inverse  to  (1-20)  only  for  the 
stated  boundary  conditions.  The  boundary  conditions  of  a  differential 
operator  should  be  considered  an  integral  part  of  the  operator,  and  :.f 
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they  are  changed,  the  operator  is'  changed.  Hence,  L**1- will  also  be 
different  if  the  boundary  conditions  of  1  are  changed. 

A  suitable  inner  product  for  electrostatic  problems  is 


01(x,y*2)  02(x,y,z)  dx  dy  dz 


(1-13) 


That  (l-13)  satisfies  (1-2),  (1-3) >  and  (l-4)  is  easily  verified.  The 
choice  of  inner  product  is  not  unique.  For  example,  the  integrand  of 
(l-13)  could  be  multiplied  by  an  arbitrary  positive  function  (weighting 
function)  w(x,  y,  z)  >0  and  it  would  still  be  an  acceptable  scalar 
product.  However,  the  particular  choice  (1-13)  makes  the  operators  L 
and  L  ^  self  adjoint,  as  we  shall  now  show. 

Let  01  and  02  represent  arbitrary  functions  in  the  domain  of  L, 
and  form  the  left-hand  side  of  (1-5) 

<$v  02>  =  rff])  $2  d-r  (1-14) 

where  dx  =  dx  dy  dz.  Green’s  identity  is 

dx  =  j^f (♦  ^  -  0  g|)  ds  (1-15) 

V  S 

where  S  is  the  surface  bounding  the  volume  V.  We  consider  S  to  be  a 

sphere  of  radius  r,  such  that  in  the  limit  r — >co  the  volime  V  includes 

all  space.  Let  0=0^  ana  y  =  0g.  The  boundary  condition  is  r0— 

as  r — **o,  which  requires  r  30/dn-**^  as  r —*■<».  Since  ds  varies  only  as 
2 

r  ,  the  surface  integral  in  ^1-15)  vanishes,  and  we  have 

fffi a^i  d*  -  JJJ 01^2  dT  (1-16) 
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Since  e  is  constant ,  (1-14)  is  symmetrical  in  0  and  0O,  and 

-L  tL 

<Upv02>  =  <0±tI0^>  (1-1?) 

We  required  02  to  satisfy  the  same  boundary  conditions  as  0^,  hence  L 
is  self  adjoint.  If  L  is  self  adjoint,  so  is  L-"1,  since 

<4v<t>2>  =  l~1p2>  (1‘l8) 

The  mathematical  concept  of  self  adjointness  of  an  operator  corresponds 
to  the  physical  principle  of  reciprocity. 

It  is  evident  from  (1-10)  and  (1-12)  that  L  and  L  ^  are  real 
operators.  We  shall  show  that  they  are  also  positive  definite,  that  is, 
they  satisfy  (1-6) .  Again  we  need  to  show  it  only  for  L  or  L”1,  since 
(1-18)  is  valid  for  =  02'  ccniu6ate  operation  is  not  needed 

since  L  is  real. )  Form 

<0,I0>  =  jJJ 0(-€  (1-19) 

and  use  the  vector  identity 

0^0  =  V  *  0V0  -  v0  *  V0  (1-20) 

and  the  divergence  theorem.  The  result  is 


e  ^  d-r 


(1-21) 


where  S  bounds  V.  Again  take  S  a  sphere  of  infinite  radius,  and  the  last 
term  of  (1-21)  vanishes  because  of  the  boundary  condition  r0-*£^  as  r — *-oc. 


<0,I0>  ^  JJje(V0.f  dt 


(1-22) 


and,  since  e  >  0,  L  is  positive  definite.  The  mathematical  concept  oi 


-5- 


positive  definiteness  is  often  related  to  the  physical  concepts  of  work 
or  energy.  In  (l-19)>  the  right-hand  side  is  proportional  to  the  electro¬ 
static  energy. 


D.  Method  of  Moments.  We  now  discuss  a  general  procedure  for  solv¬ 
ing  linear  equations ,  called  the  method  of  moments.  Consider  the  inhomo¬ 
geneous  equation 

L(f)  =  g  (1-23) 


where  L  is  a  linear  operator,  g  is  known  and  f  is  to  be  determined. 

Let  f  be  expanded  in  a  series  of  functions  fn,  ?2 ,  fy  ...  in  the  domain 
of  L,  as 


f  = 


a.f . 

i  x 


(1-2 4) 


where  the  a.  are  constants.  We  shall  call  the  f .  expansion  functions 
x  i  — - - 

or  basis  functions.  For  exact  solutions,  (1-24)  is  normally  an  infinite 
summation  and  the  f  form  a  complete  set  of  basis  functions.  For 
approximate  solutions,  (1-24)  is  a  finite  summation.  Substituting 
(1-24)  into  (1-23)  j  and  using  the  linearity  of  L,  we  have 


1 


!«(-.;)  =  g 


(1-25) 


It  is  assumed  that  a  suitable  inner  product  <f  ,g>  has  been  determined 
for  the  problem.  How  define  a  set  of  weighting  functions,  or  testing 
functions,  wn,  w„,  w,,  ...  in  the  range  of  L,  and  take  the  inner  product 
of  (1-25)  with  each  w_. .  The  result  is 

ti 

V  at.  <w_.,  Lf.>  =  <w.,g>  (1-26) 

L-.  -  3  d 

X 

j  =  1,  2,  3 ,  ...  This  set  of  equations  can  be  written  in  matrix  foxm 
as 
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where 


[V  ■ 

<wl,Lfl> 

<w1,lf2> 

*  *  * 

<w2,Lfx> 

<w2,Lf2> 

*  *  * 

*  *  o  * 

(1-28) 


_Qi 

cgj]  = 

"<Vl,g> 

°2 

<w2;g> 

m 

* 

* 

a 

(1-29) 


If  the  matrix  [i]  is  nonsingular  its  inverse  exists, 

then  given  by 


fail  =  r/-J]  [gJ 


The  Oh  are 

(1-50) 


and  the  solution  for  f  is  given  by  (1-2*4).  For  concise  expression  of  this 
result,  define  the  matrix  of  functions 


[?.]  =  fg  fj  ...]  (1-5D 

and  write 

f  =  [?,]  [a,]  =  [?.]  [iTh  tgj  (1-32) 

^  ^  -*■  tl 

This  solution  may  be  exact  or  approximate,  depending  upon  the  choice  of 

the  f .  and  w. . 
x  x 

If  the  matrix  [i]  is  of  infinite  order,  it  can  be  inverted  only  in 
special  eases,  for  example,  if  it  is  diagonal.  The  classical  eigen¬ 
function  method  leads  to  a  diagonal  matrix,  and  can  be  thought  of  as  a 
special  case  of  the  method  of  moments.  If  the  sets  f^  and  v,  are  finite, 
the  matrix  is  of  finite  order  and  can  be  inverted  by  known  methods. 


A  principal  task  in  the  solution  of  any  particular  problem  is  the 


choice  of  f.  and  w. .  She  f .  should  be  chosen  so  that  they  are  relative- 

IX  1 

ly  independent  functions  and  so  that  some  superposition  (l-2l }  can  approx¬ 
imate  f  reasonably  well.  The  v  should  be  chosen  so  that  they  are  also 
relatively  independent  and  so  that  the  products  <w .  , g>  test  relatively 
independent  properties  of  g.  We  shall  say  more  about  this  in  Section 
when  we  discuss  the  stationary  nature  of  the  solution. 


Some  additional 
(1)  the  accuracy  of 
matrix  elements ,  (3) 
(<i }  the  realization 


.  factors  which  affect  the  choice  of  f.  and  w.  are 

l  i 

solution  desired,  (2)  the  ease  of  evaluation  of  the 
the  size  of  the  matrix  that  can  be  inverted,  and 
of  a  well  conditioned  matrix  [i]. 


£.  Example ;  Charged  Conducting  Plate.-  Figure  1-1  represents  a 
square  conducting  plate,  2a  meters  on  a  side,  lying  in  the  plane  z  -  0 
with  center  at  the  origin.  Let  o(x,y)  represent  the  surface  charge 
density  on  In  ite,  assumed  infinitely  thin.  The  electrostatic  potential 
at  any  point  in  space  is 


o(x*  ,y*  5 
ajr  Tfiil - 


{1-335 


where  K  -  (x-x*)'  +  (y-y*)‘  *  -  The  boundary  condition  is  0  =  V 

(constant)  on  Use  plate.  The  Integral  equation  for  the  problem  is  there¬ 


fore 


dx’ 


izr  * 
AJ 


g(xS  y* } 


(x-x*)fc  4-  (y-y’r 


(1-3*1) 


x‘  <  a,  ;.y  a.  The  charge  density  u(x-y)  is  the  unknown  to  be  determined. 


A  parameter  of  interest  is  the  capacitance 


(1-35) 


&  *v 

Cs|-J  /"  o{x,y) 


which  can  be  calculated  once  r%x,y)  is  found. 

For  an  approximate  soluti )n,  consider  the  plate  to  he  divided  into 
H  square  subareas,  as  shown  in  Fig.  1-1.  let  the  subareas  he  denoted  * 
As,,  As r%J  . ..,  As„.  Define  functions 


on  As. 


\  =  < 


i 

V  -‘-JUy 


0  on  a  i  f  ot-ner  ii2 _ 


and  let 


j(x,y)  ~  A.  f. 


(1-37) 


where  a.  are  constants  to  he  evaluated.  In  other  words,  we  are  approxi- 

X  “  “ 

sating  a  by  a  constant  over  each  subarea.  If  (1-37)  is  substituted  into 
the  integral  equation  (1-3*0  and  the  resultant  equation  satisfied  at  the 
center  (x^  ,y_. )  of  each  As .,  we  obtain  the  matrix  of  equations 


V  -  V  i^.a. 
L  ■  ■  ir1"  ~ 


-  ij 


wnere 


'  ejix.-x'Y1  +  (y.- 


Note  that  £ . .  is  the  potential  at  the  center  or  As .  -due  to  a  unit  charge 


density  over  As . . 


solution  of  the  simultaneous  equations  (1-3-8)  gives 


an  approximation  to  the  -charge  density  on  the  plate  according  to  (1-37)  < 
The  corresponding  approximation  to  (1-35) 


t* 


H 

L  “i^i 

1=1 


gives  the  capacitance  of  the  plate. 


(i-40) 


To  pat  the  above  solution  into  the  notation  of  the  method  of 
moments,  let 


u(xjy)  *  0^x*3T/ 


(1-41) 


srf  v 


)  =  ¥ 


(1-42) 
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fCxSy*) 


4g£  v\x-x*  )2  e  (y-y')2_+  r2 


(1-43) 


Then  (1-23)  is  equivalent  to  (1-34) .  Define  the  inner  product  as 

a  a 

<f,&>  =  /*  dx  f  dy  f(x,y)  g(x,y)  (1-4%) 

-a  -i. 

The  uhkaoun  f  =  a  is  approximated  by  (l-3T)  -’here  the  ft  are  defined  by 
(1-36).  The  vei gating  f?mction3  are  defined  by 


<5,  f>  =  f 

-a 


*/ 

-a 


dy  v.(x,y)  f(x,y)  =  f(x  .,y.) 

v  J  U 


(1-45) 


These  v-  do  not  exist  as  ordinary  functions,  but  are  symbolic  functions. 

v  -  ~  _  _ 

In  particular j 

*j<x,y)  =  sCx-xJ  S(y-y^)  (1-16) 

vhere  6(x)  is  the  Birae  delta  function.  Ibv  the  elements  of  the  matrix 
(1-28)  are  given  by  (1-39),  a^d  those  of  (1-29)  =re 
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Hence#  (1-27)  is  equivalent  to  (1-38) >  and  (1-30)  represents  the  solution. 

For  numerical  results,  the  £  4  of  (1-39)  must  be  evaluated.  Let 

ji 

2b  =  2a/$f  denote  the  side  length  of  each  Aa^ .  The  potential  at  the 
center  of  As^  due  to  unit  charge  denisity  over  its  own  surface  is 


D  D 

'u-M 

-b  -b 


<Jy 


4jt€  jx  +  y2 


=  ||  (1  +  n/2)  =  ||  (0.88137) 


2b 

ire 


(l-48) 


2 

This  derivation  used  Dwight  200.01  and  731*2.  The  potential  at  the 

center  of  As .  due  to  unit  ch irge  density  over  As .  can  be  similarly 
D  1 

evaluated,  but  for  most  purposes  need  not  be.  Usually  it  is  sufficiently 
accurate  to  consider  the  charge  on  As^  as  a  point  charg . .  and  use 


,  *i 
ji  *  4jreR77 


+  (yj  -  y±)£ 


(1-49) 


This  approximation  is  3*8  percent  in  error  for  adjacent  subareas,  and 
has  even  less  error  for  nonadjacent.  ones.  Table  1-1  shows  calculations 
solving  (1-38)  and  (1-39)  for  various  numbers  of  subareas.  The  second 
column  uses  approximation  (1-49) ,  the  third  column  evaluates  £ . .  more 
precisely.  A  good  estimate  of  the  exact  capacitance  is  C/2a  =  4000 
micromicrofarads.  Figure  1-2  shows  a  plot  of  the  charge  density  along 
the  subareas  nearest  the  center  line  of  the  square  plate.  Note  that 
a  exhibits  the  well-known  square  root  singularity  at  the  edges  of  the 
plate . 

This  example  illustrates  two  simple  but  useful  approximations,  (1) 
the  use  of  expansion  functions  f.  each  of  which  exists  over  only  a 
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distance  along  plate 

Approximate  charge  density  on  subareas  closest  to 
the  centerline  of  a  square  conducting  plate. 


Table  1-1. 


Capacitance  of  a  unit  square  plate 
(micromicrofarad/meter ) 


No.  of 
subareas 

C/2a 

(approx,  .  ) 

C/2a 

(exact  f ..  ) 

1 

3150 

3150 

Q 

s 

3730 

3680 

1.6 

3820 

3770 

36 

3920 

3870 

100 

3950 

restricted  region  of  space  (method  of  subsections),  and  (2)  the  satis¬ 
faction  of  the  approximate  equation  at  specific  points  in  the  space 
(point  matching  method).  These  approximations ,  though  crude,  give 
acceptable  results  for  many  purposes. 


F.  Variational  Interpretation.  The  special  case  of  the  method  of 
moments  for  which  the  approximating  functions  f^  are  equal  to  the  weight¬ 
ing  1  "unctions  is  Known  as  Galerkin's  method.  That  Galerkin’s  method 

is  equivalent  to  the  Rayieigh-Ritz  variational  method  is  well  known. 

That  the  general  method  of  moments  is  also  a  variational  method  is  usually 
not  noted,  but  the  proof  is  essentially  the  same. 

We  first  interpret  the  method  of  moments  according  to  the  concepts 
of  linear  spaces.  Let  (Lf)  denote  the  range  of  L,  (Lf. )  denote  the 
space  spanned  by  the  Lf. »  and  (w. )  denote  the  space  .spanned  by  the 
wi-  The  method  of  moments,  Eqs.  (1-26),  equates  the  projection  of  Lf 
on  (w.)  to  the  projection  of  the  approximate  Lf  on  ,;(w, ).  Figure 
1-3  represents  this  pictorally  (f  denotes  approximation  to  f).  In  the 
special  case  of  Galerkin’s  method,  (w. )  =  ,’(f. ).  Because  the  process 


/  (LF) 


rn-iorai  representation  of  the  i 
of  moments  in  function  space. 


of  obtaining  projections  minimizes  the  error,  the  method  of  moments  is 
an  error  minimizing  procedure.  Because  the  error  term  is  orthogonal  to 
the  projections,  it  is  of  second  order.  Uhls  same  conclusion  is  obtain- 
ed  by  means  of  the  calculus  of  variations.  The  derivation  will  not 
be  given  here,  but  we  shall  summarize  the  results. 

Given  an  operator  equation  Lf  =  g,  it  is  desired  to  determine  a 


functional  of  f 

p(f)  =  <f,g’>  (1-50) 

This  functional  may  be  f  itself  if  g*  is  an  impulse  function.  Let  La 
be  the  adjoint  operator  to  L,  and  define  an  associated  function  f  *  by 


It  can  then  be  shown^  that 


_  <f,g'><r,g> 

p  <Lf,f’> 


(1-52) 


is  a  variational  formula  with  stationary  value  (1-50)  when  f  is  a 
solution  to  Lf  =  g  and  f!  a  solution  to  (1-51)-  For  an  approximate 
evaluation  of  (1-50),  let 


f »  _  V  ft  -- 

-  LP i  "i 


(1-55) 


and  substitute  into  (1-52).  It  can  then  be  shewn"*  that  the  necessary 
and  sufficient  equations  for  p  to  oe  stationary  are  those  of  (1-26), 
that  is,  the  method  of  moments.  This  variational  procedure  is  known 
as  the  Bayleigh-RItz  method,  and  hence  the  method  of  moments  is  Identi¬ 
cal  to  it. 

The  second  equation  of  (1-55)  gives  us  some  additional  insight  into 
how  to  choose  the  weighting  functions.  For  good  results  they  should 
chosen  so  that  some  linear  combination  can  closely  represent  the 
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be 


associated,  field  i 5 ,  that  is,  the  field  of  whatever  g '  appears  in  { 1-pQ ) . 
For  example,  if  we  want  f  itself,  then  g'  is  an  impulse  function,  and 
f  xs  the  Green’s  function.  This  is  a  poorly  behaved  function,  difficult 
to  approximate  oy  a  s impie  set  of  w. .  Hence,  we  should  expect  calculations 
of  the  field  itself  to  converge  less  slowly  than  calculations  of  other 
parameters  which  are  associated  with  a  well-behaved  f ' . 
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II.  FORMULATION  OF  ELECTROMAGNETIC  PROBLEMS 


A.  Introduction.  Mo  ;  engineering  problems  in  electromagnetic  theory 
consist  of  a  source,  some  matter,  ana  a  measurement  bo  be  performed. 
Antenna  system  problems  and  radar  scattering  problems  are  f vpical  examples , 
Basically,  the  solution  involves  a  determination  of  the  interaction  of 
each  element  of  matter  with  every  other  element,  with  tne  source,  and 
with  the  measuring  device.  The  equation  describing  the  problem  is  an 
inhomogeneous  operator  equation  of  the  form  of  equation  (1-1).  As  shown 
in  Chapter  I,  the  operator  equation  can  be  transformed  into  a  matrix 
eauation.  The  interaction  between  elements  of  matter  is  characterised 
..he  .a o  i.x  operator,  wi.c  Cvupj.j.n^  o _  source  lo  L-n—  matter  is 

characterized  by  an  excitation  matrix,  and  the  coupling  of  the  matter 
to  the  measuring  device  is  characterized  by  a  measurement  matrix.  An 
exact  solution  for  arbitrary  excitation  usually  requires  matrices  of 
infinite  order.  Approximate  solutions  can  be  obtained  using  matrices  of 
finite  order.  'The  solution  to  the  problem  is  represented  by  the  inverse 
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ion  cf  tns  network,  end  messuicmeiii  uf 


on  wa sponds  to  voltage  sxcitatic 
field  quantity  corresponds  to  a  weighted  combination  of  the  network  currents. 
Dielectric  bodies  are  characterised  by  two  generalized  network  matrices, 
one  independent  of  the  dielectric  permittivity  of  the  body,  and  the  other 
dependent  on  the  permittivity.  Electromagnetic  excitation  of  the  body 
corresponds  to  voltage  excitation  of  the  two  networks  connected  in  series, 
and  a  field  measurement  again  corresponds  to  a  weighted  combination  of 
the  terminal  currents.  Magnetic  bodies  are  the  dual  problem  to  dielectric 
bodies,  and  have  a  dual  representation  in  terms  of  generalized  networks. 
Hence,  they  also  are  characterized  by  two  matrices,  one  independent  of 
the  permeability  and  the  other  dependent  on  the  permeability.  Electro¬ 
magnetic  excitation  of  the  body  corresponds  to  current  excitation  of  the 
two  networks  connected  in  shunt,  and  a  field  measurement  corresponds  to 
a  weighted  combination  of  the  terminal  voltages.  Bodies  having  both  di¬ 
electric  and  magnetic  properties  are  -characterized  by  three  network  matrices, 
one  of  which  is  independent  of  the  permittivity  and  permeability,  one 
dependent  on  the  permittivity,  and  one  dependent  on  the  permeability. 

Field  excitation  of  the  body  corresponds  to  voltage  excitation  of  some 
ports  connected  in  series,  and  current  excitation  of  other  ports  connected 
in  shunt.  Field  measurement  corresponds  to  a  weighted  combination  of  the 
currents  at  those  terminals  excited  by  voltage  sources,  and  of  the  voltages 
at  those  terminals  excited  by  current  sources . 

The  classical  eigenfunction  method  of  solution  corresponds  to  obtain¬ 
ing  a  network  matrix  having  only  diagonal  elements.  The  inverse  matrix 
is  also  diagonal,  with  elements  equal  to  the  reciprocal  of  the  elements 
of  the  original  matrix-  Matrix  inversion  is  then  trivial,  and  the  infinite 
order  of  the  matrix  unimportant.  However,  the  eigenfunctions  are  known 
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If  the  network  matrix  can  be  approximated  by  a  well-conditioned  matrix  of 
finite  order,  the  inversion  is  readily  accomplished  by  known  algorithms. 

For  example,  a  high  speed  computer  can  invert  a  matrix  of  the  order  one 
hundred  by  one  hundred  in  a  few  minutes.  The  realization  of  a  well-behaved 
network  matrix  of  reasonable  order  therefore  constitutes  an  engineering 
solution  to  the  problem.  This  Chapter  deals  primarily  with  the  formulation 
of  the  desired  matrix  equations  and  their  network  interpretation. 


3-  Conducting  3odies.  Figure  2-1  represents  the  general  problem 
of  a  material  body  in  the  presence  of  an  electromagnetic  source.  For 
this  section,  it  is  assumed  that  the  body  is  a  perfect  electric  conductor. 
Let  g1  be  the  impressed  electric  field,  that  is,  the  field  produced  by  the 
source  when  the  body  is  absent.  Let  gS  be  the  scattered  field,  that  is, 
the  field  produced  by  currents  on  the  body.  The  scattered  field  is 
related  to  the  conduction  current  £  on  the  body  S  according  to 


£(*')  re(£,r’ 


)  ds’ 


(2-1) 


where  r  {r,^’ )  is  the  tensor  Green's  function  relating  a  current  element 

at  r !  to  its  electric  field  at  r.  The  total  field  at  any  point  in  space 
i  s 

is  g  +  g  and  the  boundary  condition  is  that  the  tangential  components 

,  "i  3 

£  x  (g~  g  )  =  0  on  S,  where  u  is  the  unit  normal  to  S.  Hence,  the 
appropriate  equation  for  determining  £  is 

u  x  L^(£}  =  -u  x  g1  on  S  (2-2) 

where  is  the  operator  defined  by  (2-1).  A  suitable  inner  product 
satisfying  (1-2),  (l-3),  and  (l-U)  is 
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(2-3) 


<J,E>  =  ♦  £  ds 

■r 

which  is  the  quantity  defined  as  reaction.  Note  that  (2-3)  involves 
only  the  tangential  components  of  £,  since  £  must  be  tangential  to  S. 

The  method  of  moments  is  now  applied  to  (2-2)  •  Let  the  current 
be  expanded  In  s  series  of  functions  jT^,  defined  on  S,  as 


*  T 

n 

where  the  I  are  complex  constants.  A  substitution  of  (2-4 )  Into  (2-2), 
and  an  application  of  the  linearity  of  L,  yields 


ux  Vi  L  (J  }  =  -u  x  E~ 

—  n  e'~n'  —  — 


(2-5) 


n 


are  vectors 


Nov  define  a  set  of  testing  functions  V, ,  g_,  . which 

tangential  to  S?  i.e.,  they  are  current-type  vectors.  The  method  of 
moments  requires  that  (2-5)  be  valid  for  the  Inner  product  with  each  «  x  y  , 
that  Is 


rl  <#  }L  J  >  =  -  <W  jE'S- 
n  mJ  e  n  m* 


{2-6} 


n 


for  all  m.  The  unit  normal  u  nas  been  dropped  from  (2-6)  since  and 
J  are  all  tangential  to  S- 

Ifcw  define  the  following  matrices 


Cxa3  =  j 

“  T  ~j 
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"-  <W1,SX> 
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-  <w2,e*> 
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♦ 
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( p  *?\ 
I  i 


n  3  —  I  j  ^  «L 

mn  |  X'  e  x  X'  e2 


and  rewrite  (2-6)  in  matrix  fora  as 


£z  ]  £x  j  =  [Vj 

an  n  a 


(2-9) 


The  solution  for  the  extension  coefficients  is 


vnere 


[ij  =  Cy_3  [v_  ] 


fv  1  _ 


(2-10) 


itrix  inverse  to  [2-  j.  Define  the  matrix  of  current  functions 
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elution  for  the  current  itself  can  he  »rt 
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This  equation  is  approximate  or  exact  depending  upon  the  choice  of  J ,  and 


The  designation  of  tne  sac-races  m  hne  aoove  solution  as  l: 


«Z_i,  and  iY j  is,  of  course. 
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the  corresponding  voltage,  current,  ispe-dance,  and  admittance  matrices 
of  H-nort  network  theory.  This  analogy  is  ssssrirss  hv  Fig*  2-2-  The 
conducting  body  is  characterised  by  a  generalised  in  re-dan  c  e  matrix  CZJj 

le  network  ratric 


or  its  inverse*  a  nener; 


Ci  j  - 


basically  depend,  only  on  the  geometry  or  tne 


no  tne  *3 


on  the  impressed  field  to  which  the  bodv  is  subiected.  The 
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analogy  according  to  (2-7)*  The  resultant  input  currents  of  the  network 
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[Zmn] 

depends  on  geometry 

Em 

*  and  wavelength 

• 

• 

• 

Fig*  2—2  Network  analogue  for  a  conducting 
an  arbitrary  impressed  field. 


C.  Dielectric  Bodies.  Now  let  Idle  boay  of  Fig.  2-1  be  a  dielectric 
of  Permittivity  s,  which  msy  be  ft  or  of  position  or  even  a  tensor. 

Again  the  impressed  electric  field  is  denoted  by  JJ1,  and  the  scattered 
field  by  £°.  The  scattered  field  is  produced  by  the  polarisation  currents 
J  over  the  body  according  to 

t  -  LeU)  =  Jffitz' )  re(£.£’)  dt'  (2-lU) 

V 

which  differs  from  (2-1)  only  in  that  the  integral  is  over  a  volume  V 

i  s 

instead  of  the  surface  S.  The  total  field  at  any  point  is  given  by  +  £  , 
•and  within  V  the  polarization  current  is  given  by 

si  =  jw(€  -  eo)(JBi  +  £S)  (2-15) 

s 

where  £q  is  the  permittivity  of  free  space.  Using  (2-14)  for  £  ,  one 
can  rearrange  (2-15)  into  the  form 

si 

L  (si)  -  irrr-  =  -I1  in  V  (2-16) 

ev>c/  ji»>Ae 


where  Ae  =  e  -  This  is  the  appropriate  equation  for  determining  j£. 

The  left-hand  side  of  (2-l6)  could  be  redefined  as  a  single  operation  on 
but  it  is  more  convenient  to  consider  it  as  the  sum  of  two  operations. 
A  suitable  inner  product  for  this  problem  is 


E  dt 


(2-17) 


which  is  again  a  reaction. 

For  the  method  of  moments,  let  the  polarization  current  £  be  expanded 

,  sj^,  ...,  defined  over  V,  as 


in  a  series  of  functions  sl^>  slg 


(2-18) 
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where  the  I  are  complex  constants.  Substituting  (2-18)  into  (2-16), 
and  using  the  linearity  of  L.  one  obtains 


n 

Define  a  set  of  testing  functions  I,,  over  V,  and  require 

the  equality  (2-19)  to  told  for  the  inner  products  with  each  W  .  Hence 


I 


n 


<W  ,L  J  >  -  i-  <W  ,  -2>>  =  -  <W  jES 


men 


.1“ 


m  Ae 


m' 


(2-20) 


n 


for  all  m.  Now  define  matrices  [lQ]  and  [V^j  as  in  (2-7),  and  two 
impedance  matrices,  one  [Z^3  as  in  (2-8),  and  the  other 


A 


[2mn 


3 


-1 

3« 


Ag^  .  .  .  "j 

J 

<W2,J1/fc>  <W2,J2/A£>  .  .  .  i 


(2-21) 


Equation  (2-20)  can  then  be  written  in  matrix  form  as 


[zmn  +  z  ]  [i  3  =  [v  3 

mn  mn  n  m 


(2-22) 


The  solution  for  the  expansion  coefficients  I  is  again  given  by  (2-10) 

where  ^ 

[Ym]  -  «Znm+Zwr  ]  (2-23> 

Tn  terms  of  the  matrix  (2-12)  of  expansion  functions,  the  polarization 

current  is  given  by  (2-15).  Again  the  solution  is  approximate  or  exact, 

depending  upon  the  ^  and  j^. 

In  terms  of  generalized  network  parameters,  one  can  interpret  this 
solution  as  two  impedance  matrices  connected  in  series  as  shown  in  Fig.  2-5. 
One  matrix,  [Z],  depends  only  on  the  geometry  of  the  body  and  the  wavelength. 
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Fig,  2-3  Network  analogue  for  a  dielectric  body 
arbitrary  impressed  field. 
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The  Impressed  electric  field  determines  the  voltage  excitation  of  the 

networks  connected  in  series.  The  resultant  terminal  currents  then 

determine  the  coefficients  I  of  the  oolariaation  current.  note  that  as 

n 

€ — *■<»,  the  network  analogy  of  Fig.  2-3  reduces  to  that  for  the  conduct¬ 
ing  body.  Fig.  2-2. 


P.  Magnetic  Bodies.  If  the  body  of  Fig.  2-1  is  magnetic  but  not 
dielectric,  the  problem  is  dual  to  that  for  the  dielectric  body,  discussed 
above.  To  be  specific,  let  g^  denote  the  impressed  magnetic  field,  and 

g 

g  the  scattered  field,  produced  by  the  magnetization  current  {J  over  the 
body.  The  relationship  of  gs  to  g  is  the  dual  to  (2-14)  for  the  dielectric 
case,  that  is, 


=  L  (g) 
nr*-' 


-///' 

V 


n(r) 


rm(x,£’) 


dr* 


(2-24) 


let  ji  denote  the  permeability  of  the  body,  and  u_  that  of  free  space. 
The  magnetic  current  of  magnetization  is  then  dual  to  (2-15),  or 


U  =  Mv  -  »0)(£  +  gs) 


(2-25) 


Finally,  combining  (2-24)  and  (2-25),  one  has  the  equation  for  determining 


M 


yts)  - 


in  V 


(2-26) 


which  is  dual  to  (2-l6).  Define  an  inner  product 


<M,H> 


g  dT 


(2-2?) 


The  minus  sign  is  certainly  not  necessary  in  (2-27),  but  is  used  so  that 
the  inner  product  corresponds  to  that  defined  in  the  reaction  concept. 
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To  apply  the  method  of  moments ,  let  the  magnetisation  current  g  be 
expanded  in  a  series  of  functions  gy  defined  over  V,  as 


*~T  VA 


(2-28) 


where  the  are  complex  constants, 
that  ,  „  . 


£  vn{VKn)-^rh*81 


From  (2-28)  and  (2-26),  it  follows 


(2-c9) 


.Again  define  a  set  of  testing  functions  g^,  g^,  , over  V,  and  take 

inner  products  of  (2-29)  with  each  W  .  Then,  dual  to  (2-20),  one  has 


ZV  J<W  ,L  H  >  -  4t 
n  j  m  m  n  ju 


m  Zip 


(2-30} 


for  all  m.  The  following  generalised  network  matrices  are  now  defined 


I  *  I 

L  :  J 


fV  - 


I  - 


(2-315 


[Ymn]  *  |" W  ^ W 


(2-32) 


[Y  l 

mn 


-x 


<W1^M2-/^>  *  *  ~ 

<W2,Mi/Z5p>  <Wg  ,Mp/ .  .  - 


The  matrix  equation  for  (2-30)  is  then 
[Ymn  +  Ymn]  [Vn]  = 


(2-335 


(2-34) 
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which  is  dual  to  (2-2g).  The  solution  for  the  [Y^]  i*  the  inverse  of 

(2-5^5 .  The  magnetic  current  of  the  magnetization  i*  given  by  (2-28)  once 

the  I  are  found, 
n 

The  network  representation  of  this  solution  shown  in  Fig.  2-h 

is  the  dual  case  of  Fig.  2-3.  The  matrix  [Y  1  depends  only  on  the  geometry 

mn 

A 

of  the  body  and  the  wavelength,  and  the  matrix  [Y  ]  depends  on  the 
permeability.  The  impressed  magnetic  field  determines  the  current  excitation 
of  the  networks  connected  in  shunt.  The  resultant  terminal  voltages 
then  correspond  to  the  expansion  coefficients  Vq  for  the  magnetic  current. 

B.  Bodies  both  Dielectric  and  Magnetic.  If  the  'oody  has  both  e  and 
p,  different  from  their  free-space  values,  a  combination  of  the  preceding 
two  analyses  must  be  used.  For  this  purpose,  in  addition  to  the  relation¬ 
ship  between  electric  current  and  electric  field,  (2-1  A),  and  magnetic 
current  and  magnetic  field,  (2-2^),  one  needs  the  relationship  between 
electric  current  and  magnetic  field 

E  =  Nte)  =  jfjto )  r(r,r’ )  dT*  (2-35) 

V 

A 

where  r(r,r ’ }  is  the  tensor  Green’s  function  which  relates  an  element 
of  £  at  r*  to  its  magnetic  field  at  r,  and  the  corresponding  relationship 
between  magnetic  current  and  electric  field 

£  =  -N(M)  (2-36) 

where  N  is  the  same  operator  defined  in  (2-35)*  The  minus  sign  difference 
between  (2-35)  and  (2-36)  reflects  the  minus  sign  difference  between  the 
two  curl  equations  of  Maxwell.  The  integral  equation  for  the  problem 
is  now  a  matrix  of  (2-l6)  and  (2-26),  with  the  appropriate  interaction 
terms  added.  To  be  explicit 
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Fig.  2-4  Network  analogue  for  a  magnetic  body  in  an 
arbitrary  impressed  field. 
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is  the  equation  to  be  solved  for  ;T  and  £4  vhen  the  body  is  both  dielectric 
and  magnetic. 

Equation  (2-37)  is  just  a  more  complicated  case  of  an  operator  equation. 
This  is  evident  if  the  following  matrices  are  defined 


f  = 


(2-38) 


X  = 


N 


-ii 


m- 


(2-39) 


and  (2-37)  rewritten  as 

*(f  >  +  7lM  =  -gX 


(2X0) 


Of  course,  _£  and  ?//  could  be  combined  into  a  single  operator,  and  (2-tO) 
written  as  (l-l) .  The  appropriate  inner  product  for  this  problem  is 


M*  H)  d-r 


(2 


1) 


2  ** 

which  is  precisely  the  general  definition  of  reaction.  ^ 

A  solution  by  'the  method,  of  moments  is  obtained  as  follows . 
'electric"  expansion  and  weighting  functions  as 


Define 


*n 


ana 


"magnetic"  expansion  and  weighting  functions  as 


(2-42) 


f*. 

n 


j 

r  “i 

_  !  m 

0  v  = 

— 

0  1 

j  n 

i 

-H 

-W* 

L  nJ 

n 

(2-43) 


The  expansion  for  f  is  then  of  the  fora 


f  =  V  (l  fe  +  V  f®)  (2-44) 

'  n  n  n  n'  v  3 

n 

where  the  1^  and  correspond  to  the  of  (1-24)*  Following  the  method 
of  moments,  one  obtains  the  matrix  equation 


r~n 

J 

CzL. 

i . — 

r~» 

fln] 

'fU 

f43 

+ 

— 

uh 

13 

Ec  1 

L  5211 

&m\ 

}\\ 

A*3 

(2-45) 


Here  the  various  matrices  have  the  same  definitions  as  in  sections  II-C  and 

H-B.  The  additional  matrices  fB  j  and  [C  ]  describe  the  interaction 

mn  mn 

between  electric  and  magnetic  currents,  and  the  superscript  i  has  been 
added  to  the  source  terms  (right  hand  side)  to  distinguish  them  fromt 
the  response  terms* 

The  generalized  network  representation  of  Eq.  (2-45)  is  shown  in 
Fig.  2-5.  The  network  denoted  [L  ]  again  depends  only  on  geometry 

A 

and  wavelength,  not  on  e  or  The  network  [2^3  in  series  with  the 

voltage  sources  depends  on  e,  and  the  network  [Y  ]  in  shunt  with  the 

mn 

current  sources  depends  on  u .  The  impressed  electric  field  determines 
the  voltage  sources  according  to  (2-7) ,  and  the  impressed  magnetic  field 
determines  the  current  sources  according  to  (2-31)* 


F.  Measurement.  Figure  2-6  represents  a  general  problem  of  electro¬ 
magnetic  engineering,  one  consisting  of  a  source,  a  receiver,  and  material 
bodies.  Explicitly,  Fig.  2-6  shows  a  source  connected  to  a  transmitting 
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Fig.  2-5  Network  analogue  for  o  body  having  both 
dielectric  and  magnetic  properties. 
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receiver 


Fig.  2-6  A  typical  electromagnetic  engineering  problem. 


antenna,  a  measurement  device  connected  to  a  receiving  antenna,  and  a 
scattering  object.  However,  any  of  the  material  bodies  msy  be  absent j 
herce  a  general  problem  includes  the  transmitting  antenna,  the  receiving 
antenna,  and  free«space  scattering  as  special  cases.  The  source  and 
receiver  may,  of  course,  be  at  the  same  point  in  space,  as,  for  example, 
in  an  input,  impedance  problem. 

The  voltage  at  the  receiver  terminals  due  to  a  current  source  at 

the  transmitter  terminals  can  be  expressed  in  terms  of  the  transfer 

1  3 

impedance  Z  t  between  the  two  sets  of  terminals.  '  In  terms  of  reaction, 
this  transfer  impedance  is  given  by 

Zrt  =  JJJ( Et  •  J1'  -  r  •  M1* )  dr  (2-46) 

where  E  ,H  is  the  field  produced  when  1^  is  applied  to  the  transmitter 
terminals,  and  <Tr,  Mr  is  the  current  resulting  when  Ir  is  a  lied  to  the 
receiver  terminals.  If  e  and  \i  are  scalars,  or  even  symmetric  tensors, 
the  usual  reciprocity  theorem  applies,  and  (2-46)  can  be  rewritten  as 

Zrt  =  (r  ■  -  H1  •  M*)  dT  (2-47) 

Kr,  H1  is  the  field  produced  by  1^  applied  to  the  receiver  terminals, 
and  J^,  !t^  is  the  current  produced  by  1+  applied  to  the  transmitter 
terminals.  The  important  fact  to  note  is  that  the  measurement  is  of  the 
form  of  an  inner  product  of  the  current  caused  by  an  excitation  of  the 
transmitter  with  a  field  caused  by  an  excitation  of  the  receiver.  Hence, 
in  the  general  notation  of  the  method  of  moments, 

measurement  =  <f,  g  >  (2-48) 
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1* 

where  g"  is  function  depending  on  the  measurement  to  be  performed.  For 
example,  gr  is  a  E1 ,  H1  in  (2-47).  Equation  v2«48)  is  the  general  expression 
for  a  functional  of  f .  In  terms  of  the  solution  obtained  by  the  method 
of  moments,  equation  (1-32), 

measurement  =  <f,gr>  =  [g£]  C^3  [g^j  (2-49) 


where 


r  7 


!  <fvg  >i 
r  ! 

<f2'8  >  I 


i 


J 


(2-50) 


This  measurement  matrix  is  of  the  same  form  as  the  excitation  matrix  [g^] 

of  (1-29),  except  that  the  inner  products  are  with  the  f  instead  of  the 

w  . 
n 

As  an  example  of  a  measurement  in  a  field  problem,  consider  a  conduct¬ 
ing  body  in  a  plane-wave  incident  field,  and  let  the  measurement  be  the 
tistatic  radar  echo  at  some  distant  receiver.  By  definition,  let  Echo 
be  the  field  quantity  whose  majjnitude  squared  is  the  conventional  echo 
area,  that  is 

c  =  jEchoj2  (2-51) 


A  formula  for  Echo  is  obtained  by  letting  both  the  transmitter  and  receiver 
of  Fig.  2-6  recede  to  infinity  according  to  an  appropriate  limiting 
procedure. ^  The  result  is 


Echo 


-T)  k 
2  ST  j  j 


HX 


a?  ds 


(2-52) 


where  r\  =/p/s,  k  =  2i t/k,  ^  is  a  normalized  (unit  amplitude)  plane  wave 

t 

from  the  receiver,  and  £  is  the  current  on  the  conducting  scatterer  when 
excited  by  a  normalized  plane  wave  from  the  transmitter.  The  cuirent  on 
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the  scatterer  has  already  been  determined  by  the  method  of  moments  in 


section  ii-a,  the  result  being  (2-13)* 
(2-52)  gives 


nfarrrQfTrtn  in 


Echo 
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n  nm  m 


(2-53) 


where 


3  is  given  by 


-  <<J,  ,E  > 
1-  o 


- 


(2-5*0 


This  is  the  same  form  as  the  voltage  excitation  matrix  (2-7),  except 
that  replace  the  In  the  Galsrkin  method  =  j^)  the  two  voltage 

matrices  [v^J  and  [Vnl  are  just  different  excitations  of  the  same  generalized 
network.  In  the  method  of  moments,  [i7]  may  be  interpreted  as  the  excitation 
of  an  adjoint  generalized  network,  defined  as  the  transpose  of  [z^]* 

G.  Discussion.  The  method  of  moments  is  so  general  that  almost  any 
solution  can  be  interpreted  as  an  application  of  it.  The  use  of  eigen¬ 
functions  for  the  expansion  functions  and  adjoint  eigenfunctions  for  the 
weighting  functions  is  the  classical  eigenfunction  method.  Point-matching 
procedures,  that  is,  satisfaction  of  the  operator  equation  at  specific 
points  on  a  body,  arc  equivalent  to  using  impulse  functions  as  weighting 
functions  in  the  method  of  .nomer  s.  Methods  which  approximate  the  operator, 
for  example,  the  method  of  nets,  are  interpretable  as  special  cases  of 
the  general  moLiiod.  Hence,  one  should  view  the  method  of  moments  as  a 
unifying  concept  rather  than  as  a  particular  technique  of  solution.  In 
applying  the  method  to  engineering  problems,  one  is  confronted  with  many 


possible  choices  and  approximations.  The  particular  technique  of  approxi¬ 
mation  called  the  method  of  subsections,  introduced  in  Section  I-E,  has 
been  found  quite  useful  for  digital  computation. 

Basically,  the  method  of  subsections  can  be  thought  of  as  dividing 
the  object  into  a  number  of  pieces,  and  calculating  the  interactions 
according  to  the  method  of  moments.  One  of  the  first  applications  to  time- 


harmonic  electromagnetic  fields  was  the  calculation  of  scattering  by 

h 

square  conducting  cylinders.  General  computer  programs  for  basically 
the  same  method,  but  with  somewhat  better  approximations,  have  been  used 

for  cylinders  of  arbitrary  shape.  Solutions  of  this  type  have  also  been 

6  7 

obtained  for  dielectric  cylinders,  for  conducting  bodies  of  revolution,1 

g 

and  for  wire  antennas  of  arbitrary  shape.  An  analysis  of  wire  antennas 
and  scatterers  of  arbitrary  shape,  both  loaded  and  unloaded,  are  formulated 
using  the  methods  of  this  section  in  the  next  chapter. 

In  terms  of  the  method  of  moments,  the  method  of  subsections  involves 
using  expansion  functions  which  each  exist  over  separate  sections  of  the 


body,  and  using  weighting  functions  which  test  the  field  in  each  section 
(often  by  point  matching).  In  the  case  of  conducting  bodies,  this  is 
equivalent  to  obtaining  the  self-impedance  of  each  subsection  of  current, 
and  the  mutual  impedance  of  each  element  with  every  other  element  of  current. 
In  the  case  of  dielectric  (or  magnetic)  bodies,  the  interpretation  is 
similar,  except  that  the  self  impedance  (or  admittance  in  the  magnetic 
case)  can  be  divided  into  two  components,  one  dependent  on  e  (or  u)  and 
the  other  independent  of  €  (or  p) .  The  mutual  terms  are  always  Independent 
of  e  (or  p).  Hence,  the  careted  networks  (left-hand  ones)  of  Figures  2-3, 
2-4,  and  2-5  become  noninueracting  elements,  that  is,  expressible  as  diagonal 
matrices.  For  example,  the  network  representation  of  Fig.  2-5  becomes 


that  of  Fig.  2-7  when  the  method  of  subsections  is  used-  The  effect  of  e 


independent 
of  €  and  fjL 


and  \i  is  then  just  a  loading  of  the  appropriate  terminals  of  the  general¬ 
ized  network  for  the  body. 

The  use  of  these  techniques  for  computation  gives  rise  to  a  number 

of  practical  difficulties,  which  can  be  overcome  by  careful  analysis. 

For  example,  the  field  of  a  discontinuous  distribution  of  current  has 

singularities,  and  the  order  of  differentiation  and  integration  cannot 

always  be  interchanged.  The  usual  tensor  Green’s  function  often  cannot 

be  used,  but  must  be  interpreted  as  symbolic  of  the  more  fundamental 

o 

vector  and  scalar  potential  formulas.  The  accuracy  of  a  solution  and 
the  ease  of  computation  depend  to  a  large  extent  on  the  ingenuity  and 
care  used  in  the  formulation  of  the  problem.  The  method  of  subsections 
is  basically  a  small  body  technique,  since  the  matrices  involved  become 
too  big  for  large  bodies.  The  treatment  of  large  bodies  requires  great 
care  in  the  choice  of  expansion  functions  and  testing  functions,  and 
even  then  the  interaction  matrix  usually  is  valid  for  only  a  subclass 
of  excitations.  Ferhaps  the  best  approach  to  the  large  bodies  is  a  per¬ 
turbation  one,  that  is,  the  use  of  an  approximate  solution  as  one  of  the 
expansion  functions.  The  method  of  moments  then  gives  a  correction  to 
this  approximate  solution.  A  disadvantage  to  this  approach  is  that  the 
generality  of  the  solution  is  lost,  being  valid  for  only  one  excitation 
of  the  body.  Much  more  work  remains  to  be  done  on  the  large  body  problem. 
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III.  WIRE  .ANTENNAS  AND  SCATxSRERS  OP  ARBITRARY  SHAPE 


A.  Introduction.  A  general  analysis  ol*  objects  constructed  of  thin 
conducting  wires  ar.d  situated  in  an  arbitrary  impressed  electromagnetic 
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analysing  thin  wire  antennas,  -as  discussed  in  Chapter  II.  However, 
for  any  particular  problem  there  are  many  different  ways  to  obtain  a 
solution.  The  conventional  retarded  potential  formulation  is  chosen  here 
because  of  its  conceptual  simplicity-  The  resultant-  equation  differs 


r.HP  yyg 


iX> 


SO  lilt 


T*r 


■ _ iivi  i_A_  lii  l 


msHH  :  - 


the  wire  axis.  (J)  The  boundary  condition  (3-5)  is  applied  to  the  axia. 
component  of  jg  at  the  wire  surface-  To  this  approximation ,  Eos.  (3-1)  to 
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wi*n  e auations  similar  to 


The  s*  s  are  given  is  terns  of  the  I’s 


(5-12)  and  (5-15)  lor  g(s)  and  o(5}. 

by  (5-15).  and  hence  (5-10) 
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represented  by  (5-10)  as  the  equations  for  an  H  port-  nei*jcrk  with  terminal 
pairs  (a,  a) .  Be  voltage  applied  to  each  port  is  approximately  |** 

Hence,  by  defining  matrices 


M  l  = 


!  *  i 
i  i{h)  i 


r  ? 


IV  j  =. 

I 

{ 

i 


«  i 


l(2)  -  4i 


(3-lM 


j^CS) 


5 


(3*351 


v  substituting  (5-11 


='  in's 


one  can 


Into  (5-10)  and  rearranging  into  the  fo^a  of  (5-1;?)-  Alternatively, 
aaslv  1 3-10 1  to  13-13)  to  two-  isolated  eiisests  ana  obtain  the 


n  n  05  OS$u  vcC-SSSS 


is  ssssvhat  easier  to  folio* 


a  o  oivMjyi 


consider  two  representative  elements  of  the  vise  sCc,u  i/Crcr  j 


W4*W  •*** 


in  Pig.  3-2.  The  integrals  in  (3-11)  and  (3-12)  are  of  the  same  form,  and 
are  denoted  by 


-jkR 


mn 


*(n,m)  -  -gr  j  TZ 

n  .  „ 


A* 


- di 

R  n 
mn 


(3-16) 


n 


Symbols  +  and  -  are  used  over  m  and  n  when  appropriate.  Evaluation  of  the 
\jr  in  general  is  considered  in  Section  III-F.  Let  element  n  of  Fig.  3-2 
consist  of  a  current  filament  l(n),  and  two  charge  filaments  of  net 
charge 

q(n)  =  l(n)  q(n )  =  ^  l(n)  (3-17) 


where  q  =  ga£-  The  vector  potential  at  m  due  to  l(n)  is,  by  (3-11), 


A  =  dl(n)  £&n  v(n,m)  (3-13) 

The  scalar  potentials  at  in  and  in  due  to  the  charges  (3-17)  are,  by  (3-12), 
0(m)  =  JL  [l(n)  *(n,m)  -  I(n)  iK'n.m)] 

(3-19) 

0(m)  =  ^  [I(n)  )|f(n,m)  -  l(n)  ^(n,in)] 

Substituting  from  (3-18)  and  (3-19)  into  (3-10),  and  forming  Z  = 

mn 

]<r(m)  *  A*m/l(n),  one  obtains 

Zmn  =  t(n,m)  +  ^  [>J-(n,in)  -  t-(n,in) 

-^(n,m)  +  ^(n,m) }  (3-20) 


This  result  applies  for  self  impedances  (m  =  n)  as  well  as  for  mutual 
impedances . 

The  wire  object  is  completely  characterized  by  its  impedance  matrix, 
subject,  of  course  to  the  approximations  involved.  The  object  is  defined  b 
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"°  Segmems  Of  0  Wire  severer 


2B+1  points  on  the  wire  axis,  plus  the  wire  radius .  The  impedance 
elements  are  calculated  by  (3-20),  and  the  voltage  matrix  is  determined 
by  the  impressed  field,  according  to  (3-1^)-  The  current  at  N  points  on 
che  scatterer  is  then  given  by  the  current  matrix,  obtained  from  the 
inversion  of  (3-15)  as 

[I]  =  [Yj  [V]  [Yl  =  [Z]-1  (3-21) 

Once  the  current  distribution  is  known,  parameters  of  interest  such  as 
field  patterns,  input  impedances,  echo  areas,  etc.,  can  be  calculated 
by  numerically  evaluating  the  conventional  formulas . 


C.  Wire  Antennas.  A  wire  antenna  is  obtained  when  the  wire  is  excited 


by  a  voltage  source  at  one  or  more  points  along  its  length.  Hence,  for 
an  antenna  excited  in  the  n-th  interval,  the  applied  voltage  matrix  (3-3>) 


is 


[Vs]  =  -o 


(5-22) 


i.e. ,  all  elements  zero  except  the  n-th.  which  is  equal  to  the  source 
voltage.  The  current  distribution  is  given  by  (3-21),  which  for  the  [V] 
of  (3-22)  becomes 


(I]  =  V. 


n 


*1 


in 


2  n 


Nn 


(5-23) 


Hence,  the  n-th  column  of  the  admittance  matrix  is  the  current  distribution 
for  a  unit  voltage  source  applied  in  the  n-th  interval.  Inversion  of 
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the  impedance  matrix  therefore  gives  simultaneously  the  current  distribu¬ 
tions  when  the  antenna  is  excited  in  any  arbitrary  interval  along  its 

length.  The  diagonal  elements  Y  of  the  admittance  matrix  are  the  input- 

on 

admittances  of  the  wire  object  fed  in  the.  n-th  interval,  and  the  Y  „  are 

mn 

*he  transfer  admit : -snces  between  a  port  in  the  n-th  interval  and  one  in 
the  n-th  interval. 

The  radiation  pattern  of  a  wire  antenna  is  obtained  by  treating  the 
antenna  as  an  array  of  N  current  elements  l(n)£E^.  By  standard  formulas, 
the  far-zone  vector  potei  cial  is  given  by 


V  a 

A  =  y  I(^)  e 

r\  *—• 


KX  COS  g_ 

n  vi 


(3-21) 


where  r  and  r  are  the  radius  vectors  to  the  distant  field  point  and 
~r>  ~n 

to  the  source  points,  reflectively,  and  g  is  the  angle  between 
Xq  and  r_.  The  far- zone  field  components  are 


E_  =  -iw; 


E0  - 


(3-25) 


where  Q  and  0  are  the  conventional  spherical  coordinate  angles. 

An  alternative  forisula  for  the  radiation  pattern,  more  convenient 
for  computation,  can  be  obtained  by  reciprocity.  Figure  3-3  represents 
a  distant  current  element  I£^  (subscripts  r  denote  "receiver”),  adjusted 
to  product  the  unit  plane  wave 


E  =  u  e 


— -s»c  *  r 

n 


(3-2c) 


in  the  vicinity  of  the  antenna.  Here  u  is  a  unit  vector  specifying  the 

v  -T* 

polarization  of  the  wave,  k  is  a  wave  number  vector  pointing  in  the 


direction  of  travel  of  the  wave,  and  r  is  the  radius  vector  to  a  point 
n  on  the  antenna.  By  reciprocity,1 


I  &iJ 


(3-27) 


antenna 


where  is  the  33^  component  of  £  from  the  antenna,  and  I  is  the  current 
on  the  antenna.  The  constant  2/1*  is  that  needed  to.  produce  a  plane 
wave  of  unit  amplitude  at  the  origin,  which  is 


-j*~r 


to  u  e 


(3-28) 


i  numerical  approximation  to  (1-2  7)  is  obtained  by  defining  a  voltage 


Vf]  =  jf(l)  •  ££,  | 

—  i 


I  r(2)  *  a t2 ; 


(3-29} 


£r(N)  *  Al-r  j 


where  is  given  by  (3-26),  and  expressing  (3-2?)  as  the 


e  matrix  product 


-.jxr 


(riiTivf: 


(3'30) 


where  [V  j  denotes  the  transpose  of  [Vj. 

Note  that  hr  ]  is  the  same  matrix  as  (3' 1^5*  that  is,  it  is  the  voltage 
matrix  for  plane-wave  excitation  of  the  wire.  Equation  (3-30)  remains 
valid  for  an  arbitrary  excitation  (V  j;  it  is  not  restricted  to  the 
single  source  excitation  (3-22). 

The  power  gain  pattern  for  the  u  component  of  the  radiation  field 


is  given  oy 


?(M3  =  - 


4«ro 


(3-31) 


where  p  =  .We  is  the  intrinsic  impedance  of  space,  and  P.  is  the  cower 
*  *■  *  -  an 

input  to  the  antenns  (*  denotes  conjugate; 


=  Be  [[Vs]  [I*]}  =  Re  {[Vs]  [Y*j  [VS*]} 


0-32} 


For  the  special  c^se  (3-22)  of  a  single  rource,  becomes  simply 
Re(jVnj2  Y  }.  Using  (3-30)  and  (3-32)  is  (3-31)*  one  has 


g(9}0) .  3&  itrst&LmjnL 

43  Be  {[Vs]  [P-J  [Vs*]} 


(3-33) 


where  £v~(8,0}3  is  given  by  (3-2§)  for  various  angles  of  incidence  9,0. 
If  the  total  power  gain  pattern  is  desired,  the  g's  for  two  orthogonal 
polarisations  may  be  added  together. 


D.  Wire  Seatterers.  Consider  now  the  field  scattered  by  a  wire 
object  in  a  plane  wave  incident  field.  Figure  3-L  represents  a  sc&tterer 
and  two  -distant  current  elements.  It.  at  the  transmitting  point  £«.,  and 

i*  V 

11  at  the  receiving  point  The  li  is  adjusted  to  produce  a  unit 

oiane  wave  at  the  seatterer 


U  ~  a 


v  J 


where  the  rotation  is  analogous  to  that  of  ($-26).  The  voltage  excitation 
matrix  (3—1^)  is  then 


Cv*3  =  i  JE^Cl)  -  4£._  ! 


i  is  i 2)  *  i 
1  "2  | 


(3-55) 


,  Ay  I 

<  V-’)  A6;  I 


and  the  current  [Ij  is  given  by  (3-21)  with  [Vj  •=  [V'*].  The  field  produced 
by  [I]  can  then  be  found  by  conventional  techniques . 

The  distant  scattered  field  can  also  be  evaluated  by  reciprocity,  the 


same  as  la  the  antenna  case.  A  dicole  Jl  at  the  receiving  point,  is 


The  solution  is  in  the  form  of  an  impedance  or  admittance  representation 
of  the  scatterer.  Hence,  the  general  procedures  for  determining  scatter- 

O 

ing  from  loaded  bodies,  available  in  the  literature,  can  be  used  directly. 
This  means  that  the  behavior  of  wire  antennas  and  scatterers  loaded  by 
lumped  impedance  elements  in  any  of  the  N  intervals  can  be  calculated 
from  the  present  solution  by  matrix  manipulations.  Furthermore,  the  effect 
of  finite  conductivity  of  the  wire,  or,  more  generally,  arbitrary  boundary 
conditions  at  tne  wire  surface,  can  be  analyzed  by  considering  the  wire 
to  be  appropriately  loaded  in  all  intervals  along  its  length. 

Fv  Evaluation  of  An  accurate  evaluation  of  the  scalar  \jr  function 
of  (3-3.6)  is  desired.  Let  the  coordinate  origin  be  located  at  the  point 
m,  and  the  path  of  integration  lie  along  the  z  axis.  Then 


where 


t(m,n) 


dz' 


2a  = 

n 


(3-39) 

(3-4o) 


and  a  =  wire  radius.  The  geometry  for  these  formulas  is  given  in  Fig.  3-5. 

One  approximation  to  the  fs  can  be  obtained  by  expanding  the 
exponential  in  a  Maclaurin  series,  giving 


a 

l  mn 


(3-42) 


The  first  term  is  identical  with  the  static  potential  of  a  filament  of 

charge.  The  second  term  is  independent  of  R  .  Hence,  a  two  term 

mn 
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Fig.  3-5  Geometry  for  evaluating 
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approximation  to  (3-39)  is 


♦(“>») 


Bjta 


log 


fs  +  a  +  Jr 


+  (z+a)2 


z  -  a 


+  /p2  + 


(z-a  Y 


dk 


(3-43) 


If  r 


-/ 


2  2 

p  +  z  is  large,  then 
e"dkr 


( 


r3-44) 


For  a  first  order  solution,  one  can  take  (3-43)  as  applying  for  small  r, 
say  r  <  2a/3,  and  (3-44)  for  large  r,  say  r  >  2a/3* 

For  higher  order  approximations,  more  rapid  convergence  can  he  obtained 
by  taking  a  phase  term  e  out  of  the  integrand.  Then 


,-.1kr  r  .  w  '”mn 


8  -  r) 

/  - 
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Term  by  term  integration  gives 

-jKr 


f(m,n)  = 


ojtq: 


[I,  -  jk(l  .  n  ) 


mn 


+  ... )dz' 


(3-45) 


p 

-  —  (I  -  2rT  4-  »-  T  ) 
2  vx3  2  ‘  ~L' 


+  3  r-  (14  -  5rI3  +  5r% 


r\)  + 


J  (3-46) 


-'here 


I,  =  log 


IT 


z  +  a  +  +  (z  -f  a)* 
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/2  * 
p  +  (z  -  a)c 


(3-47) 
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(3-^8) 


Qtf-z  /  2  ,  ,  » 

-p-  vP  +  (QH-Z! 


I2  +  £?  +  (z~af 


T  X1 


(>49) 


T  2  ,  2cP  +  002 

I.  =  20b  +  - * - 

**  '  ^ 


(3-5C) 


An  expansion  of  the  type  (3-46)  is  theoretically  valid  for  all  r,  hut  it 
fails  numerically  for  large  r  because  it  involves  subtractions  of 
numbers  almost  equal  in  magnitude.  For  p  <  a,  one  should  set  p  =  a  in 
the  expansion. 

An  expression  suitable  for  large  r  is  obtained  by  expanding  (3-39) 
in  a  Maclaurin  series  in  2'  as 


~  /  [f(0)  +  f  ’  (0)z*  +  fw(0)  (zl)2  +  ...  j 


(3-51) 


where 


f(z’)  =  - 


P 


2  +  (z  -  z'}2 


p  *  (z  -  z’)2 


When  a  five  term  expansion  of  (3-51)  is  integrated  term  by  term,  there 


results 


-jkr  _  2 
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[Ac  +  jkaA1  +  (ka)  Ag  +  j(ka)^  +  (ka)  AL]  (3-b2) 


wnere 
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5<i>  f-1  +  ^  +  15  (|)5  f3 


30(j)2  +  35{p)43 


*2  -  -  5^>2  '  ^(f)2  t1  -  “(|)2  + 
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(3-53) 


\  120 


For  accuracy  of  better  than  one  percent,  one  can  use  (3-46)  for  r  <  xOO. 
and  (3-52)  for  r  ->  10a. 

An  alternative  derivation  of  the  type  of  (3-52)  can  oe  obtained  as 
follovs.  For  r  >  z‘.  one  has  the  expansion 
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where  are  the  spherical  Bessel  ^unctions  of  the  first  kind,  tr”'(kr) 
are  the  spherical  Henkel  functions  of  the  second  kind,  and  P^(z/r)  are 
the  Legendre  polynomials.  If  (3-5*0  is  substituted  into  (3-3v)  and 
integrated  term  by  term,  there  results 


t(m,n)  = 
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Equation  (3-55)  can  be  rearranged  into  the  form  of  (3-52)  -  although  the 

(2'J 

recurrence  formulas  for  rr  1  and  P  make  commutation  directly  from  (3-5t=) 

y*|  -n  *  v  \  ^  i 


axmost  as  easy. 
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TV.  CALCULATIONS  FOR  LINEAR  WIRES 


A.  Plecewlse-iinear  Current,  Galerkin’s  Method.  A  linear  wire  antenna 
or  scatterer  is  one  constructed  of  a  straight  piece  of  wire,  or  of  several 
colinear  pieces  of  wire.  Again  the  wire  cross  section  is  assumed  circular, 
and  Its  diameter  small  compared  to  wavelength.  This  is  a  special,  case 
Gf  the  problem  considered  in  Chapter  III,  and  some  calculations  wore  made 
using  that,  theory.  However,  to  obtain  faster  convergence,  the  ft ml  calcu¬ 
lations  re  made  using  a  piecewise- linear  approximation  to  the  current. 

This  is  equivalent  to  using  "triangle1’  functions  as  expansion  functions 
in  the  method  of  moments .  For  this  section,  the  weighting  functions  were 
taken  to  be  the  same  triangle  functions,  and  hence  the  solution  corresponds 
co  Galerkin's  method. 

Figure  1  shows  a  straight  section  of  wire,  and  defines  the  coordinate 
system.  The  wire  extends  from  z  =  0  to  z  =  L  along  the  z  axis,  and  is 
of  radius  a.  It  is  assumed  that  only  the  axial  component  of  the  current 
on  the  wire  is  significant,  and  it  is  expressed  in  terms  of  the  net  current 
l(z)  at  any  po.int  z.  As  discussed  in  Chapter  HI,  the  problem  is  represented 
by  the  operator  equation 


Ll(z)  =E1(  z) 


(1-1) 


where  IT  is  the  z  component  of  the  Impressed  electric  field  at  the  wire 
surface.  The  operator  L  is  determined  by  the  usual  vector  potential  method. 
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The  boundary  conditions  for  the  current  are 


1(0}  =  I(L)  =  0 


(k-h) 


since  no  cnarse  can  accumu 1  art 


Lne  enas  oi  tne  wire. 


“ox*  convenience,  consider  the  wire  to  be  (n-rl)  meters  long.  Let 


I(z)  be  apprc-xime 


;ea  ov  a  series  o 


>f  n  triangle  functions 


I(s)  =  )  I.  T(z  -  i) 


?  1  -  jzj ,  m  <  1 


T(a)  =  :: 
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;  o,  jz)  >  i 

\ 

Note  that  (t-5)  satisfies  the  boundary  conditions  (--*•).  The  approximate 
current  (t-5)  has  the  value  l(i)  =  L,  and  is  linear  on  each  interval 
defined  by  successive  integers.  Following  Galerkin’s  method  (the  method 
of  moments  with  w.  =  f . } ,  one  obtains  the  usual  matrix  equation 


>  Z..  I4  =  V. 
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Comparison  of  the  input  impedance  of  o  center  -  fed  iineor  onteono,  L/2a  »  74.2 
Present  theory _ ,  XingS  solution  ooo,  Hu's  solutions  ♦♦♦  . 


Figures  4-5  to  4-7  illustrate  how  the  current  distribution  on  the 


antenna  changes  es  the  feed  point  changes.  Only  the  magnitude  is  plotted, 
but  the  real  and  imaginary  parts  were  also  calculated.  In  each  case  the 
source  strength  was  one  volt.  Figure  4-5  is  for  I  =  X,  Figure  4-6  is 
for  L  =  l.pX,  and  Figure  4-7  is  for  L  =  24  The  case  1  =  X/2  showed 
little  change  in  the  shape  of  the  current  distribution  as  the  source 
was  moved,  and  hence  was  not  plotted. 

Figures  4-8  tc  4-jjO  show  how  the  gain  pattern  for  the  linear  antenna 
changes  as  the  feed  point  changes.  These  curves  were  normalized  so  that 
they  represent  power  gain  over  an  isotropic  radiator.  Figure  4-6  is  for 
L  =  X,  Figure  4-9  is  for  L  =  l.pX,  and  Figure  ^-10  is  for  L  =  2X.  The 
case  L  =  X/2  was  not  plotted  since  the  gain  pattern  is  relatively 
insensitive  to  feed  position.  Note  how  the  antenna  tends  to  behave  as 
a  traveling  wave  anaemia • as  the  feed  is  moved  to  one  end.  This  is 
particularly  evident  in  the  L  =  2X  case. 


D.  Linear  Wire  Scatterers.  Computations  for  linear  wire  scatterers 
excited  by  plane  waves  at  various  angles  of  incidence  have  also  been  made 
for  L/2a  ratios  varying  from  10  to  2000,  and  L/X  values  up  to  2.1.  The 


following  results  for  the  case  L/2a  =  .2  (u  =  IO)  are  giv¬ 


en  nere  to 


illustrate  the  general  behavior  of  wire  scatterers. 

Figure  4-n  shows  echo  area/X1-  vs.  L/X  and  compares  the  result  with 

2 

the  solution  of  Y.  Y.  Hu.  professor  Hu's  solution  is  actually  an  impedance 
type  solution  using  Galerkin’s  method,  with  the  wire  divided  into  two 
sections.  On  each  section  a  constant  plus  a  sinusoid  was  used  for  expansion 
functions.  As  can  be  seen  from  Fig.  4-11,  the  accuracy  of  her  solution 
is  good  ever  the  range  of  L/X  considered.  This  indicates  that  long  wires 


probably  can  be  ac._  .rately  treated  if  divided,  into  segments  of  length  up 
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Current  (mogr.itude)  on  a  linear  ontenna,  L/Za  «  74.?„  as  the  feed  point  vwies 


os  the  feed  point  varies. 


ontennc ,  •**?  *74.2 ,  L  “  I.5X  os  rh«  fwtd  point  vortos. 


<r/X 


to  a  wavelength,  if  constants  plus  sinusoids  are  used  as  expansion  functions - 
However,  such  a  procedure  gives  rather  difficult  integrals  to  evaluate, 
and  the  solution  of  this  caper  is  easier  to  compute  if  the  wire  is  not  too 


long. 

Figures  4-12  to  t-14  _jov  how  the  current  distribution  on  the  wire 

changes  as  the  angle  of  incidence  of  the  plane  wave  excitation  changes.  For 

these  curves  the  magnitude  of  the  component  of  S  parallel  to  the  wire  was 

one  vrlt  per  wavelength.  Only  the  magnitude  of  the  current  is  shown,  but 

the  real  and  imaginary  components  were  also  calculated.  Figure  4-12  is 

o  --O 

for  the  case  L  =  4,  with  angles  of  incidence  varying  from  /}  to  -O 
in  steps  of  lp°.  Figure  4-15  is  the  corresponding  set  for  L  =  1-5X,  and 
Figure  4-.lt  is  for  L  =  24.  The  current  on  a  wire  L  -  4/2  was  relatively 
insensitive  to  angle  of  incidence,  and  is  not  shown. 

Figures  4-ip  to  4-1'/  show  the  bistatic  radar  cross  section  patterns 
for  the  same  cases  as  the  current  was  shown  ir  Figures  4-12  to  4-14,  The 
angle  of  incidence  is  shown  in  each  case  by  an  arrow.  Hots  that  there 
is  a  large  lobe  at  an  angle  of  scatter  equal  to  the  angl^  of  incidence, 
out  on  the  other  side  of  the  direction  normal  to  the  wire.  This  corres¬ 
ponds  to  specular  diffraction  from  the  wire,  and  is  more  pronounced  as 
the  vi e  becomes  longer.  Figure  4-15  is  for  wires  of  length  1-  =  4,  Figure 
4-l6  for  =  1.54,  and  Figure  4- 17  for  L  =  24.  Again  the  case  L  =  4/2 
is  not  shown  because  the  shape  of  the  cross  section  patterns  is  relatively 
insensitive  to  the  angle  of  incidence. 


£.  Loaded  Wire  Scatterers.  The  impedance  elements  characterizing 
a  linear  wire  are  basically  the  same  parameters  as  the  impedance  elements 
used  in  the  theory  of  loaded  scatterers.^  Hence,  it  is  relatively  easy 


to  calculate  the  behavior  of  a  wire  with  lumped  impedances  at  points  alont: 


3L* 


e)  Angle  of  incidence  =  30* 


f }  Angle  of  g>cidincs*l5* 


Fig.  4-12  Carrot  (rnogrufuJe)  on  a  linear  cnfenrtG  ^  =  74.2  ,L4. 

for  plane  »j«  excitation (IE,I = one  vort  per  wavelength } 
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a)  Angle  of  incidence  =  90°  b)  Angie  of  incidence  *75* 


t/»  L/2  3L«  U  L/4  L/Z  3L* 


e)  Angie  of  incidence  =  30°  f)  Angie  of  incidence  =  15° 

Fig. 4-13  Current  ( magnitude)  on  a  linear  antenna -=-=74.2  ,L  =  I.5X, 

*  .  .  2a  ’  * 

for  plane  wave  excitation  (IEzl*one  volt  per  wavelength) 
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a)  Angle  of  incidence  =  90° 


b)  Angle  of  incidence  =  75 


c)  Angle  of  incidence  =  60° 


e)  Angie  of  incidence  =30° 


d)  Angle  o€  incidence  =  45 


f)  Angle  of  incidence  =  15° 


Fig. 4-14  Current  (magnitude)  on  o  linear  antenna^- =  74.2 ,  L  =2 X, 
for  plane  wave  excitation  ( IEZ I  =  one  volt  per  wavelength) 


*)Ar<gfc  of  ifxyderc* 5  30* 


Fig.  4- 15  SWalic  rodor  cross  s*ction  poft 


f  )kr*$t  of  inoc*nc**l5' 


on»»ol-^«74.2tL*X,o»«»  ongis  of  nexfenu  var*s. 


o)  Angie  of  haidence  >90* 


WAn*e  of  K*J* nce«73* 


•) Angie  of  ineiJene#  *30*  f) Angie  of  Wfwe »  «* 


R3.4-I6  B«)etic  rodor  croee  teeHon  pcHem*  for  0  Sneer  arternc,-^j-»74.2,;.»I.S\,o*  file  ongle  of  incidence 


vt 


Rg.4-17  Bsfotic  radar  cross  s*ciion  pc tlws  *oro  iintar  a«*rvxj,-^- *74.2 , t*21;cs  the  crxj»  of  nbd*x*  rones 
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its  length.  Again  calculatio  ns  far  '/2a  ratios  of  from  K)  to  2000,  and 
for  L/X  ratios  up  to  2.1,  were  made.  The  following  results  for  the  case 
1/2  a  =  7^-2  (C  =  10}  are  representative  of  the  general  behavior  of  loaded 
scattorers. 

p 

Figur-r  4-lS  shows  the  echo  area/X-  for  a  center  loaded  dipole  with 
various  loads  ZT ,  chosen  to  correspond  to  those  used  by  Professor  Hu 

jj 

2 

for  her  calculations.  Again  professor  Hu’s  results  are  good,  but 
apparently  become  worse  as  — *•<».  However,  this  discrepancy/  is  not  as 

bad  as  it  seems,  since  the  echo  area  of  an  open-circuited  dipole  is  very 
sensitive  to  the  capacitance  across  the  gap.  For  an  infinitesimal  gap, 
the  capacitance  is  infinite,  .just  as  it  is  in  linear  antenna  theory.. 

In  any  approximate  solution,  the  gap  capacitance  depends  on  the  approxi¬ 
mations  made.  In  Professor  Hu's  solution,  no  function  capable  of 
representing  a  singularity  in  charge  distribution  is  used,  hence  her  result 
give  a  low  gap  capacitance.  In  the  method  of  subsections  used  for  the 
calculations,  the  gap  capacitance  increases  with  the  number  of  subsections 
chosen.  This  is  because  the  function  expansion  (--5)  can  come  closer 
to  representing  a  singularity  as  the  number  of  subsections  increases . 
Actually,  a  very  small  adjustment  of  the  gap  capacitance  in  Professor 
Hu’s  results  would  bring  them  into  close  agreement  with  the  results  of 

this  report,  even  for  the  case  2_ — >«. 

Jj 

Figure  4-19  shows  the  echo  area/X  for  a  center  loaded  dipole 
resonated  by  reactances  at  various  values  c-f  L/X.  The  general  theory 
of  resonant  scatterers  is  available  in  the  literature.  By  definition, 
a  resonant  loaded  scatterer  is  one  for  which  the  load  impedance  is  a 
reactance  equal  to  the  negative  of  the  input  reactance  to  the  scatters.* 


when  fed  as  an  antenna.  This  definition  gives  a  maximum  echo  area  when 


the  open-circuit  echo  from  the  scatterer  is  much 


smaller  than  the  short- 


circuit  echo-  Such  a  condition  is  usually  met  for  small  scatterers, 
but  not  for  ones  of  dimensions  X/2  or  greater.  Figure  4-19  illustrates 
this  in  the  following  way.  Curve  (a)  is  for  a  center-loaded  da pole 
continually  tuned  to  resonance.  Curves  (fc)  and  (c)  are  for  the  dipole 
tuned  to  resonance  by  an  inductor  at  L/X  =  0-25  and  0.35 >  respectively. 
Curves  (e)  and  (f }  are  for  the  dipole  tuned  to  resonance  by  a  capacitor 


at  L/X  =  0.55  and  0 . 69 -  The  fact  that  curves  (e)  and  ( f }  overshoot 
curve  (a)  indicates  that  resonance  defined  in  terms  of  input  impedance 
does  not  necessarily  give  maximum  echo  area.  Curve  \d)  is  for  an  open- 
circuited  dipole.  The  curve  for  a  short-circuited  dipole  would  lie  ?.n 
between  curves  (c)  and  (e). 
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V.  SCATTERING  BY  CONDUCTING  CYLINDERS 
OF  ARBITRARY  SHAPE 


and  is  computed  using  the  large  argument  formula  for  the  Hankel  function 


Since  p  >  >  p'  we  approximate 


'  &'|  =  P  ~  P*  cos  (0  -  0’ } 


where  0  and  0*  are  defined  in  Fig-  5-1.  Then 


Le(^  =  ^  !  j  }  e 


jkp*  cos  (0  - 


r  / 


i 

ua.S 


! contour 


for  an  incident  electric  field  of  strength  tp  The  scattering  width  may 
also  he  calculated  using  reciprocity  as  illustrated  in  Section  III-D  foi 
wires  of  arbitrary  saaoe. 


3-  False  Function  Approximation  to  Current.  simplest  solution 

to  (5-1)  is  one  in  which  the  integral  over  the  surface  is  approximated 
as  a  sum  of  integrals  over  H  small  intervals.  In  each  interval  the 
current  is  assumed  uniform  ana  the  field  is  matched  at  one  noint  within 


eacn  interval. 

Figure  5-1  illustrates  the  division  of  a  symmetrical  body  into  ri 
intervals  and  defines  the  notation.  The  i-th  segment  is  given  by  the 
coordinates  of  ^ . .  at  which  point  the  integral  equation  is  satisfied. 

The  starting  and  end  no-inis  c-f  the  interval  are  defined  by  the  coordisat 
of  p.  and  p. .  The  length  of  the  interval  is  defined  as 


=  a  J-.  -I-  A  £, 


(5-M 


snssre  A  f.  —  *  ^ 


and  a"  £. .  = 


In  terms  oi  the  language  oi  tne 
functions  are  pulses  and  the  weight!: 


SU I 


metnoa  oi  moments  the  expansio: 


_G5_ 

s*— 


(5-?) 


L^C')  -  £j£(f>) 


^  2- 


(5-6) 


within 


0  outside  a£ . 

i 


r.(P.)  =  J  -  A’j)  dks. 


07) 


The  resulting  set  of  equations,  assuming  an  incident  plane  wave  of 
strength  r,  from  x  =  -oo  is 


f  a.  - 

L.  ji  i 


j  =  1,  2,  N 


(5-8) 


where 


v-y  "Oh  -“ii5  «*! 


(5-9) 


-  jkx . 


(5-10) 


The  diagonal  matrix  elements  must  be  integrated  analytically 
since  the  .ntegrand  is  singular.  This  integration  is  illustrated  in 
Section  V-E.l.  When  the  field  and  source  points  are  different  (j  i) 
no  r-mgule  cities  are  involved  and  to  a  crude  approximation 


.94- 


(5-11) 


=  h:  Mkln.  -  «  I }  k^£. 
Ji  o  ■  j  ~x'  X 


These  are  the  approximations  it  is  believed  K.  K.  Mei  made  in  his  solution. 
For  symmetric  bodies  with  normal  incidence 


ai  =  °N+l-i  a  ~  %+i-i  ^i  ^N+l-i 


(5-12) 


so  that  the  matrix  equation  (5-8)  may  be  reduced  to 


N/2 

V  a.  =  he 

L  j1  1 


5  =  1,  2,  N/2 


(5-15) 


where 


i  +  i 
ii  i,N+l-i 


=  *ii  +  Ho  ■ '<*k  ‘  pR+l-il)  k^i 


(5-14) 


ii.  *  i..  +  i.  . 

jx  ji  j,N+l-i 


-  -  %l> +  Ho2)(kk  -  **5 


(5-15) 


Note  there  are  the  same  number  of  points  on  both  halves  of  the  scatterer 
so  that  K  is  even. 

Inversion  of  the  matrix  equations  (5-8)  or  (5-15)  gives  the  cr  which 
represent  the  current  density  within  the  i-th  interval.  This  current 
may  also  be  found  by  computing  the  magnetic  field  outside  the  conductor. 


Hence 


■mv  ■  sx  *i x  hi *  +  vj  x  (5-i6) 


where  n.  is  the  normal  to  the  conductor  at  the  j-th  point.  This  method 


requires  calculating  two  matrices,  the  column 


matrix 


whose  elements  are 
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,  -jkx . 

H1  *  v  £.  x  ^  x  2  e 
.1  K  .1  0 


g 

and  the  square  matrix  IL .  whose  elements  are 

J1 


(5-17) 


!ji  =  5§k  x  Vj  x  2Ho2)<k  I  -  %l  > 


(5-l8) 


where  z  is  the  unit  vector  in  the  z  direction.  The  matrix  eauation  for 


is  then 


fJj  -  [in  +  [hs]  [a] 


(5-19) 


The  singular  matrix  elements  HT.  are  evaluated  in  Section  V-E.2,  using 
an  analytic  integration.  Equation  (5-18)  is  in  itself  a  very  crude 
approximation.  However,  it  is  easily  evaluated  using  analytic  differentia¬ 
tion  of  the  argument  since  no  singularities  are  involved.  The  rectangular 
components  of  the  normal  are  easily  approximated  as 


nx  -  *<»!  - 

n  -  (t  -  t  )M 


(5-20) 


The  solution  was  first  programmed  using  (5-^3)  and  (_;-!!)  for  the 
matrix  elements.  The  results  for  the  current  density  and  echo  area  are 
shown  in  Figures  5-2a  and  5-21 .  Note  that  even  with  such  crude  approxi¬ 
mations,  the  echo  width  as  calculated  using  the  a.'s  is  in  good  agreement 
with  Andlreasen's.  The  current  density  as  calculated  by  the  a. ' s  is 

considerably  in  error.  However,  as  calculated  by  n  x  jj,  J  is  in  good 

s 

agreement,  being  in  error  only  in  the  region  of  unequal  intervals  and 
rapid  curvature  of  the  ellipse.  No  more  current  calculations  were  made 
using  this  method.  The  reason  for  the  more  accurate  current  evaluation 
is  probably  related  to  the  fact  that  this  procedure  is  somewhat  of  an 
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iteiation  procsss  where  only  an  approximation  to  L  has  been  made. 

For  a  better  evaluation  of  the  2^.  (j  i)  of  (5-9) >  the  addition 

0- 

theorem  for  Hankel  functions'* 


(k!^  -  &’l ) 


=  £  Vm(kt)  42)(^|)  cos  m  (5-21) 


€  =1  m  =  0 


2  m  4*  0 


is  used  to  expand  the  kernel  of  the  integral  equation  around  each  in 
terms  of  the  variable  of  integration  t’  along  the  scatterer  contour. 
Figure  5-3  defines  the  notation  where  0*  is  now  the  angle  between  t '  and 


a  ~  *  Gi' 


Each  interval  is  taken  to  be  flat  so  that 


cos  0’  = 


&  *  a 


i'  ■  %  *  % 


(5-22) 


The  result  of  the  integration  for  the  2 . .  is  given  in  Section  V-E.3-  An 
upper  limit  of  four  seemed  sufficient  for  convergence  of  the  integrated 
series  (5-^9)* 

It  should  be  noted  that  the  elements  and  2.  .  are  most  efficiently 

computed  at  the  same  time  since  they  involve  the  same  values  of 
/o\ 

^  '(kl^..  -£-!)•  The  symmetry  of  the  scatterer  may  also  be  taken  into 

m  0  i 

account  using  (5-12)  and  (3-15) • 

When  the  scatterer  surface  curves  rapidly,  the  matrix  element  (5-^9) 
may  be  altered  by  defining  angles  0’  and  0’  as  shown  :n  Fig.  5-^-  The 
diagonal  elements  remain  as  in  (5-^3)  but  for  j  /  i,  the  matrix  elements 


become 
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Fig.  5-4  Illustration  defining  angles  and  j>  on 
scatterer. 


e  kA  Z.  in+if  ,  -«  KA  ii-i  /0x 

V =  2Zt"1)'°  sr  t-r^)  I  s-tt  -  ^rrT57T<-T^)  co5  ^ 


kA  £.  2' 


m=G 

4 


c-*  e  kA  i.  eh-1  f  , 

+2I  s?  Hr±> 

a=0 


®  -r  1  (m+i) 


,  kA  £  eL  T  /  ^  + 

H5557  (-5^)  I  4  Wcl)  cos 


where 


cos  0’  = 


S.’  ’  ii 


a.  -  iu 


£U 

X 


(5-23) 


>24) 


'Vfc! 

4. 

and  a  similar  equation  for  cos  0 ’ . 

Figure  5-5  shows  the  solution  for  the  a. ‘ s  for  three  different  cases. 

In  all  three,  (5-43)  is  used  as  the  singular  point  element .  The  solution 

given  oy  dots  used  (5--9)  for  the  off  diagonal  eluents  as  did  the  solution 

given  by  circles.  In  xhe  latter,  more  points  were  added  about  s  =  .5  where 

the  ellipse  curves  rapidly.  The  solution  given  by  crosses  used  (5*23) 

for  off  diagonal  elements.  In  all  three  solutions,  a  “jump"  occurs 

where  the  spacing  between  points  changes  drastically.  This  problem  remained 

throughout  the  work  and  is  probably  due  to  an  inadequate  integration 

procedure.  Furtner  work  will  be  done  to  eliminate  this  problem.  Echo 

area  has  not  been  shown  since  the  results  plotted  as  in  Fig.  5-2b- 

A  further  calculation  was  done  where  (5-43)  was  used  foi  tne 

singular  matrix  elements  and  (5-11)  for  the  £..  when  kj^j  >  1.3*  For 

J  *■ 

the  £..  where  kj&j  <  1-5?  the  better  approximation  (5-23)  was  used.  The 
31 

results  plot  identically  with  Fig.  5-^>  illustrating  that  the  current 
density  at  a  point  is  highly  dependent  on  the  field  near  that  point  and 
not  very  sensitive  to  the  field  far  away. 
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Fig.  5-5  Pulse  approximation  to  current  density  on  ellipse. 


ai+i 


Fig.  5-6  Pictorial  Representation  of  Linear  Constrain!  on  af . 


The  e.l  o  ents  of  the  reduced  matrix  may  be  calculated  as 


'  T5 
2  " 
mn 


j  C.U 


rA  ~ 

-&  g.1  i-a 


2n-2 


£  + 
2m-l,2n-l 


'll  <en 


2  a 


m  =  1,  . . . ,  M  n  =  2, 


M-l 


i,R- 

mn 


V 


2m- 1,1 


n  =  1  m  =  1,  . . . ,  M 


•R  _  o  /k-l,(N/2)-lAk(N/2)- 

»M  2m-l,(N/2)  ^ 

^(N/2)- 


n  =  M  m  =  1,  — ,  M 


(5-28) 


and  those  of  fE1]  are 
z 


i  “^^x2nv-i 
E1-  ne  d  “ 
z 


m  =  1,  . . 


M 


(5-29) 


for  a  plane  wave  from  x  =  -oo. 

The  result  of  using  (5-^.5)  and  (5-22)  for  the  f...  matrix  elements 

V 

and  (5-28)  and  (5-25)  ^or  the  reduced  matrices  is  shown  in  Fig.  5 -7. 

The  field  has  now  only  been  matched  at  19  points  but  the  result  are 
quite  good 

It  is  also  possible  to  expand  the  current  density  in  a  Taylor  series 
about  each  in  terms  of  the  body  surface  coordinate  t.  Keeping  only 
the  linear  term  in  t 


sls(ft)  =  3X3  (&)  + 


~3l - 


a 


t  = 


a.  +- 
1 


itl 


-ai-l 


(5-50) 


where  the  derivative  term  is  approximated  as  shown.  Integration  over 
one  interval  then  gives  a  contribution  to  the  i  -  1,  i,  and  i  +  1  elements 
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3lO. 
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Fig.  5 -7  Linearly  constrained  current  density  an  ellipse. 


of  a  row  of  the  matrix.  A  matrix  element  is  not  computed  directly 
but  is  obta?  .led  as  a  result  of  integrating  over  three  intervals  on  the 
scatterer  surface  and  summing  the  proper  contributions.  The  contributions 
to  the  matrix  elements  of  the  j-th  row  as  a  result  of  integrating  over 
the  i-th  interval  (column)  is 


i  /  Ho2' Wj-aJP  +  f  H^Oelaj-oil)  dkt: 

(5-31) 


*i 


where  £&.  has  been  defined  in  (5-25).  The  first  integral  has  been  computed 

-i. 

previously  for  both  j  =  i  and  j  £  i  and  the  results  given  by  (5-43)  and 
(5-23).  The  second  integral  is  evaluated  in  Section  V-E.4. 

A  piecewise  linear  approximation  to  the  current  density  may  also 
be  made.  In  this  case,  triangles  are  chosen  as  expansion  functions  in 
the  method  of  moments  and  the  weighting  functions  are  again  delta  functions 
Thus  N 

is  -  Y.  QlTi  (5-32) 


1=1 


wnere 


1  + 


-A  k .  ^  t  <  0 


",  =  < 
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Integration  over  the  i-th  interval  then  contributes  to  the  matrix  elements 
of  the  j-th  row  as 
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(3-34) 
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The  integrals  here  have  all  been  computed,  when  doing  a  Taylor  series 


approx.Lmui..Lo:i  to  Current  density.  ±n  tins  ease,  noWever,  tne  matrix 

elements  have  been  constructed  from  somewhat  different  contributions. 

For  symmetric  scatterers  use  is  again  made  of  (5-15) • 

The  dots  in  Fig.  5-3  give  the  result  for  the  Taylor  series  approxi¬ 
mation  and  the  circles  are  the  result  for  the  piecewise  linear  approximation 


D.  Solution  Using  Galerkin's  Method.  A  solution  for  the  TM  polariza¬ 
tion  was  also  attempted  using  the  expansion  and  weighting  functions  to  be 
the  same  in  the  method  of  moments. 

Cnoosing  i  pulse  approximation  to  the  current 


w.  =  f.  =  P.  =  < 

i  2.  I  1 


within 


0  outside  AS. 


(5-55) 
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where 


2..  =  /  P^(t")  LP..  ( t * )  dkt! 


(5-57) 


with  LP. (t1)  given  by  (5-7)-  Here  primed  variables  denote  "source" 
coordinates  and  double  primed  variables  denote  field  coordinates.  For 
integration  purposes  the  Hank  el  function  with  j  j  i  is  expanded  in  both 
"source”  and  "field”  variables  in  Section  v-il-5-  The  interval  is 
considered  "flat"  as  in  fig.  5-5-  In  the  resulting  expansion,  integration 
of  the  small  argument  Bessel  function  yields  terms  identical  to  those 
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Fig.  5-8  Tayior  series  and  piecewise  linear  approximations  to  current  density 
on  ellipse. 


Puteo  approximation  to 
method . 
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and  the  H&nkei  functions  are  evaluated  at  x’  =  ,  y *  =  y .  ,  x  =  x., 

—  tJ 

y  =  y . .  Tills  is  exactly  a  Taylor  series  in  four  variables  including 
*  r* 

second  derivative  terms.  The  integrations  involve  powers  of  (x  -  x. ) 
and  (y  -  y. }  and  are  performed  in  a  fashion  similar  to  that  used  to 
evaluate  (5*38). 

The  result  of  programming  this  method  plots  very  close  to  the  solution 


given  in  Fig.  P-9*  There  are  several  shortcomings  to  this  method.  The 


expand cd  equations  contain  many  terms,  although  most  are  easily  programmed. 
However,  a  piecewise  linear  approximation  to  the  current  density  requires 
terms  in  vxsth  variables  x  and  y.  Inspection  of  Galerkin’s  method  indicates 


•.hat  the  resulting  number  of  terms  in  the  integration  would  be  very 
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Evaluation  of  Integral  s . 
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[i.. j  with  Poise  Approximation  to  Current .  The 


integration  is  over  a  sman.  interval  so  that  tne  smau.  argument 
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may  be  used.  Now  t 1  is  the  variable  of  integration  along  the  i-tn  segment 


of  the  scatterer  with  origin  at  the  point  defined  o\ 


Hence  (p-9) 


becomes 


k£  i. 
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2.  Evaluation  of  fji?  J  for  Poise  Approximation  to  Current.  7c  evaluate 
s 

B. .  a  small  argument  approximation  to  the  Hankel  function  say  again  be 
used.  However  since  the  singularity  is  nonintegrabie,  the  integral  is 
first  expressed  with  the  field  point  off  the  surface  of  the  conductor 
a  small  distance  n.  Thus 
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On  integration 
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4.  Evaluation  of  Linear  Term  j.n  Expansion  for  Current.  The  integration 
is  done  exactly  as  in  Sections  V-E'.l  and  V-E-3*  For  the  singular  terms 
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For  integration  over  nonsingular  intervals 
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function  is  cow  written  as 


(i  -  t^n  icg  2%Ht+  (i^-)2 


for  t  >  t ’  and  with  t  and  t  *  interchanged  for  t  <  t  * .  The  integration 
to  he  carried  out  is  then 
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After  the  first  integration.’  is  carried  out  the  second  may  he  reduced 
considerably.  The  final  result  is 

iu  -  Ci  -  (t^i)2) 

i  (k A£.)2  (3-2  log  (r  — 2~)(1  '  ^i)2)  “  m  ^i>2) 

(5-57) 


F.  References 


1.  K.  K.  Mei  and  J.  Van  Bladel,  "Scattering  from  Perfectly  Conduct¬ 
ing  Rectangular  Cylinders,"  IEEE  Trans . ,  vol.  AP-11,  March,  1963, 

pp.  185-192. 

2.  M.  G.  Andreasen,  "Scattering  from  Parallel  Metallic  Cylinders 
with  Arbitrary  Cross  Sections,"  TREK  Trans.,  vol.  AP-12,  November,  19&4, 
pp.  7^6-75^- 

3.  R.  F.  Harrington,  "On  the  Calculation  of  Scattering  by  Conducting 
Cylinders,"  IEEE  Trans.,  vol.  AP-13>  September,  i$65,  pp-  812-813- 

4.  R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  McGraw- 
Hill  Book  Co.,  New  York,  19&1,  p.  23<?. 

5.  J.  A.  Stratton,  Electromagnetic  Theory,  McGraw-Hill  Book  Co., 

New  York,  1$&1 ,  p.  37^ • 


-118- 


VI.  SCATTERING  3Y  DIELECTRIC  CYLINDERS 


A.  Introduction.  In  this  section,  the  problem  of  scattering  by  di¬ 
electric  cylinders  is  discussed  and  a  particular  solution  is  formulated 
in  terms  of  electric  and  magnetic  surface  currents  on  the  obstacle.  The 
fields  inside  arid  outside  are  represented  in  terms  of  the  field  radiated 
by  electric  and  magnetic  surface  currents  distributed  over  the  surface 
of  the  obstacle. 

The  two  dimensional  curve  representing  the  surface  of  the  obstacle 
is  approximated  with  straight  line  segments  and  currents  are  assumed  to 
be  constant  over  each  straight  line  segment  or  to  be  delta  functions 
located  at  the  center  of  the  strip.  The  fields  due  to  infinite  strips 

or  ^  directed  either  in  the 
2  direction  or  along  the  straight  line  segment  are  calculated.  A  dual 
formulation  is  used  to  extend  the  results  to  the  TE  polarization. 


oi  electric  or  magnetic  current,  with  ^ 


B.  Formulation  of  the  Problem  (IM  to  z).  Consider  the  problem 
indicated  in  Fig.  6-1.  A  plane  wave  with  electric  field  in  the  z  direction 
(TM  to  z)  is  incident  upon  an  infinite  dielectric  cylinder  of  arbitrary 
cross  section,  peirittivity  e,  and  permeability  jjl. 

We  can  represent  the  fields  in  terms  of  incident  and  scattered  fields 


E  =  ES  +  Ein° 
z  z  z 


w3? 


me 


(6-1) 


H  =  F*  +  HinC 
P  'p  P 

For  purposes  of  calculating  the  fields  external  to  the  cylinder  we 
can  use  the  equivalence  principle,  replacing  the  induced  polarization 
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Fig.  6-1  Cytinderical  S cofferer  ond  Coordinate  System 
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ana  magnetisation  sources  b£/dt  and  v  x  £1  with  surface  currents  ,T  =  q  x 


and  =  £  x  |i  distributed  over  the  surface  of  the  obstacle.  These  sources 
are  such  as  to  produce  zero  fields  inside  (if  the  currents  are  exact). 

t 

Since  these  equivalent  sources  produce  zero  fields  inside,  we  can  replace 
tt  cylinder  vivn  *?*/  material  and  the  external  fields  will  not  be  affected. 
If  we  choose  the  material  as  free  space,  then  the  problem  becomes  a  homo¬ 
geneous  one  of  calculating  the  fields  due  to  ^  in  homogeneous 

free  space. 

Similarly,  for  purposes  of  calculating  the  internal  fields,  we  can 
replace  the  effects  of  the  incident  field  with  surface  currents  -J  s= 

-Hi  x  £  and  =  *£  x  n  distriouted  over  the  surface  of  the  cylinder. 

These  sources  are  such  as  to  produce  zero  fields  outside  the  cylinder 
(if  the  currents  are  exact),  so  that  we  can  replace  the  free  space  outside 
the  cylinder  with  any  material  and  the  internal  fields  will  not  be  affected. 
If  we  choose  the  material  to  be  of  permittivity  e  and  permeability  p,  then 
the  problem  becomes  the  homogeneous  one  of  calculating  the  fields  due 
to  -H  in  infinite  homogeneous  space  of  characteristics  p,e. 

Vector  potential  functions  can  be  calculated  for  these  two  homogeneous 
problems .  Denoting  potentials  outside  by  subscripts  1,  then 


A,  =  ri  S  JU’)  H^2'(k 


J 

contour 


is(^t)  Ho  (k0P}  * 


r  C. 


=  ij  9  -t2)(v) 


contour 


a*’ 


(6-3) 


P  =  /(x  -  x‘)  +  (y  -  yT 

df  *  =  differential  element  of  length  around  contour  of 
cylinder 

k  =  m  ,/e  u 
o  v  o*o 
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'i'he  fields  outside  are  given  by 


fl.  =  a1  +  v  x  &  +  -A—  V  x  V  x  F 
1  1  *1 


(6-h) 


The  vector  potentials  for  the  fields  inside  are  similar  except  the  material 
wave  j.umher  and  material  constants  armear.  Tins 


*s  -  SJ  /-  HoS)(kop) 


3  (£')  H^(ke)  &£' 
“s'  o  v 
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22  =  -Vx£2  +  ._S?xVx&! 

%  =  v  x  +  kt  v  *  v  x  e2 


where  k  =  ^Jue 

Since  there  are  no  surface  currents  in  the  actual  problem*  the 
tangential  jg  and  g  fields  must  i  e  continuous  at  the  surface  of  the 


cylinder.  On  the  surface  of  the  cylinder 

a  x  a  =  a  x  % 


(6-6) 


Equations  (£-6)  may  be  expressed  using  (6-4)  and  (6-5)  in  terms  of  two 
integro-differential  operator  equations 


a  x  i- 


Vx,7x4  +  iS?VxVx4  +  Vx£i*Vx£2J=ax^ 


a  X  [-7  X  ^  ♦  v  X  4  -  ji-  v  X  V  x  El  ♦  4j  7  *  V  X  Eg]  =  a  x  rf* 


(6-7) 
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which  are  of  the  form; 


Ln(fi)  +  -  *1 

^21^1^  +  i22‘I2;  =  % 


‘  /  Q\ 

kc-0) 


i>-.  Lp,  •  Lpp  are  integro-differential  operators,  f  1  f  are 


the  unknowns  ^  and  g,  ,gD  are  the  known  functions  n  x  ?Tnc  and 


a  >'  a 


Equations  (6-7)  say  be  rewritten  in  the  fora 


r_j  _j  jn  j 

a  4‘4J  =  a x £ 

r.JT  .  j  .J4  Ml  i 

ax[V%  +  %'%J=^xfe 


me 


(6-9) 


i 

where  J3p  ic  the  field  inside  du_  to  the-J^  source  in  infinite  homo geneou s 
media  with  constants  p  and  e  and  g,  is  the  field  outside  due  tc  the  J 

x  *-S 

source  in  free  space.  Like  definitions  hold  for  the  g  fields. 

In  order  to  convert  the  operator  equations  to  a  set  of  matrix 
equations  the  contour  of  the  scatterer  is  divided  into  N  straight  line 
segments  as  in  Fig.  6-2.  The  current  is  then  assumed  uniform  within 


eac a  interval  so 


N 


*• -  I  “i  pi 

i=l 

3 

l-lhh 


(6-10) 


vnere 


wiiGin  i— th  segment 


?.  = 


I  0  outside  i-th  segment 

i 

Using  (c-iO)  and  satisfying  (6-9)  at  the  center  of  each  segment  gives 
two  sets  of  matrix  equations.  For  3M  polarization 
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Magnetic  currents  care  assumed  counterclockwise . 
Electric  currents  are  assumed  in  1  direction . 


Fig. 6-2  Dielectric  scctterer  and  notation . 
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To  illustrate  the  form  of  the  four  sou  are  submatrices 
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where  each  matrix  element  is 


■Jj  ~J  _  ^ 


0X  <*1 


(6-1-) 


Here  the  partial  matrix  element  Sg..  is  the  electric  field  at  the  center 

of  the  j-th  segment  due  to  a  current  of  magnitude  unity  in  the  i-th 

interval  when  the  medium  has  permeability  y  --ad  permittivity  e.  is 

■*■*/■* 

the  same  type  of  field  in  free  space.  Similar  definitions  hold  for  the 
other  three  submatrices  and  their  matrix  elements.  A  dual  set  of  equations 
holds  for  the  TE  polarisation  as  discussed  in  YT-3-  'These  are  obtained  by 
substituting  the  dual  quantities  in  (6-11)  to  (6-14).  The  remainder  of 
this  chapter  is  then  devoted  to  finding  the  matrix  elements  for  iM 
ix>  lari  cation . 


C.  Evaluation  of  Partial  Matrix  Elements  From  Magnetic  Currents,  we 
wish  to  calculate  the  fields  at  the  midpoint  of  the  j-th  segment  due  to  the 
unit  magnetic  current  at  the  i-th  segment  of  Fig.  6—2. 
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It  is  assumed  that  the  magnetic  currents  are  all  directed  counter 


clockwise*  Then,  locating  ourselves  at  the  midpoint  of  a  segment  and 
facing  in  the  direction  of  the  magnetic  current  arrow,  the  fields  for 
the  subscript  1  case  (free  space)  will  be  just  to  the  right  and  the 
fields  for  the  subscript  2  case  will  be  just  to  the  left. 

Consider  the  case  i-j ,  that  is  the  fields  are  to  be  evaluated  at  the 
midpoint  of  the  j-th  segment  due  to  magnetic  currents  in  the  j-th  segment. 
The  orientation  of  the  coordinate  system  is  unimportant  for  this  problem 


and  the  configuration  may  be  simplif'  ed  to  that  of  Fig.  6-5- 

The  magnetic  field  is  obtained  from  the  vector  potential  ]?,  vher< 


* 3 


E^2'(ko)  dx* 
o' 


(6-15) 


-d/2 


and 


Sow  is  dielectric  media  with  constants  4  and  e 


fi  = 


V  x  V  x  Z 


( 6- 16) 


Since  in  Fig.  6-5  the  tangential  magnetic  field  and  the  partial 


matrix  element  Hi  -  _  are  identical,  then 


*(? 

j  E^(kp)  dx’ 

J  o 

^  -d/2 


(6-17) 


This  field  may  be  evaluated  in  the  following  way.  let  the  field  point 
lie  on  the  y  axis,  at  a  point  (0,y),  evaluate  /  IT*-'  (kp)  dx’,  differentiate 
twice  with  respect  to  y,  and  then  let  y — wO.  Hence 


Fig.  6-3  Magnetic  Current  Source 
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h!L  =  H  (0,0+)  =  Ha 

y — >0 


sd 
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(kp)  ax'] 


Taking  the  first  tvo  terms  of  each  series  in  the  series  expansion  of 

(2^  2 
Hv  "(kp),  using  integration  formulas  Dwight  623  and  Dwight  623*2, 
c 

differentiating  and  taking  the  limit  as  y— ►£,  the  partial  matrix  element 
becomes 


3  t  7kd 

2  "  log 


(6-19) 


To  calculate  the  tangential  magnetic  field  H  (0,0)  in  the  subscript 

A 

1  case,  consider  the  same  limiting  process  with  the  field  point  approaching 
(0,0)  along  the  negative  y  axis;  by  symmetry  the  tangential  &  field 

will  be  identical  to  that  obtained  in  the  subscript  2  case,  since  the 
magnetic  current  source  is  the  negative  of  the  source  in  the  subscript 

2  case.  Thus  we  obtain  the  field  by  replacing  the  unit  magnetic 
current  with  its  negative  and  p.,  e,  k  with  eQ,  kQ.  Hence 


ri?. .  =  -  J-  (k  a)  -  +  — 


3  ,  *oa 

5- 106  -r 


(6-20) 


The  tangential  electric  field  on  both  sides  of  a  strip  of  magnetic 
current  of  strength  unity  must  also  be  obtained.  From  Fig.  6-3,  if  we 
assume  that  the  electric  field  to  the  left  of  the  midpoint  is  in  the  2 
direction,  then  from  symmetry  considerations,  the  electric  field  just 
to  the  right  is  in  the  negative  z  direction.  The  strength  of  the  field 
is  most  easily  obtained  by  using  the  dual  of  Ampere’s  law 
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(6-21) 


d£  =  K 


where  K  is  the  magnetic  current  enclosed.  The  result  for  a  unit  magnetic 
current  is  easily  obtained  as 


20J 


(6-22) 


Now  consider  the  evaluation  of  the  fields  when  j  0  i.  Replace  the 

magnetic  surface  current  distribution  with  a  magnetic  dipole  element  of 

strength  d.  located  at  the  midpoint  of  the  segment.  Here  d,  is  the 
1  •• 

length  of  the  i-th  segment.  For  purposes  of  calculation,  assume  a  magnetic 
current  element  located  at  the  origin  and  oriented  in  the  x  direction  as 
in  Fig.  6-3-  Then  to  a  crude  approximation 


£  -  sj  h£2)(kp) 


E  =  -V  x  £ 


Thus 


2  =  4-VxVx£ 


kd, 

*3 

kd 


kd .  /  p  \  i 

=  FT"  £in  0 


H  -_i  SSll  h(2)'  (ho) 
p  4q  kp  o 


(6-23) 


(6-24) 


kd. 


I U  =  - 


sin  0  (kp) 


Now  if  we  locate  the  magnetic  current  element  at  the  center  of  the 
i-th  strip  with  an  orientation  corresponding  to  that  of  the  i-th  strip 
we  can  derive  fields  at  the  center  of  the  j-th  strip  by  replacing  p 


with  |P.-|  and  0  with  0...  -  0. .  Then 
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(6-25) 


v 

kL  =  — rrF  sin  (04J 

JA 


-  0j  H$?i)(k|p,,i ) 


J.  ■  Ji- 


4i  =  ‘V  Si”  (S*J  '  V  +  Hi=  005  C0J  ‘ 


-kd. 

x 

¥tT 


(2) "/ 


sm^  -^.)  sin  -  0,)  '  (k|Pai|) 


cos  (5#  -  0  )  cos  (0  -  0.) 

- TTp^ - iL^-«i)NPjii) 

10 


(6-26) 


The  elements  E^f. .  and  E?. .  can  be  obtained  by  replacing  k,  e,  p. 
xj  x  x  3 1 

in  (6-26)  with  k  ,  e  ,  u  . 
v  '  0  0  *0 


D.  Evaluation  of  partial  Matrix  Elements  from  Electric  Currents.  We 
wish  to  calculate  the  fields  at  the  j-th  point  due  to  a  unit  electric 
current  at  the  i-th  segment. 

First  consider  the  case  when  j  =  i  as  in  Fig.  6-3.  Then 


/V 

J  o 


(2) 


(kp)  dx' 


so 


d/2 


=  /  ^2)(kx')  dx' 


(6-27) 


(6-26) 


Since  E  and  E„ . .  are  identical,  we  have  on  integrating  directly 
z 


_  -kc  i  d  (2),kdi  nd  (2)  ,kd,  ,kd,  „(2),kd.  ,kdv  l 
-  33£  ]  2  Ho  +  T  [H1  (_2}  VT?  Ho  H1^}J  ? 

1  (6-29) 


(?) 

where  H.  and  H  are  Struve  functions .  Alternately,  if  Hv"'(kx')  is 

X  o  o 

approximated  by  its  series  expansion  (retaining  the  first  two  terms  of 
each  series),  and  integrated  directly,  the  following  result  is  obtained 
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\ 

) 


E233  '  -  (“X1  '  3<“) 


(6-30) 


The  element  is  obtained  by  replacing  K  by  Kq  in  (6.50). 

The  tangential  magnetic  field  may  be  obtained  by  applying  Ampere’s  law 


in  a  manner  similar  to  that  used  in  Section  VJ-C  to  calculate  the 
tangential  electric  field.  The  result  is 


=  ICL, 


(6-3 1) 


Now  consider  the  fields  when  j  /  i.  Replace  the  electric  surface 
current  ^  with  an  electric  filament  directed  in  the  2  direction  with 
current  I  =  a. ,  and  located  at  the  midpoint  of  segment  i.  Consider  a 
filament  located  at  the  origin.  Then 


A  =  Jr  JT2'(kP) 

a  ii  -■  a  \  * 


Hence 


r(2)  i 


E  =  — jJL  ;(ko) 

2  o  v  7 


(6-32) 


(6-53) 


Hp  =  0 

The  partial  matrix  elements  for  the  subscript  2  case  are  then 


(k| 


piii) 


Sill 


K 


(6-3*0 
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The  elements  and  are  obtained  by  replacing  k  and  rj  by  kQ 
aid  in  (6-34) .  The  matrix  elements  in  (6-11)  have  thus  essentially- 
been  established. 

E.  Duality.  For  incident  fields  TE  to  z,  the  problem  may  be  solved 
by  duality.  In  this  case,  all  magnetic  fields  will  be  in  the  z  direction 
and  the  electric  fields  will  lie  in  the  xy  plane .  Thus  the  scattered 
fields  can  be  calculated  by  assuming  magnetic  current  sheets  in  the  z 
direction  and  electric  currents  sheets  directed  along  the  contour  of 

the  scatterer  in  the  xy  plane.  Again  the  problem  is  resolved  into  two 
homogeneous  problems,  one  involving  free  space  and  the  other  involving 
an  infinite  homogeneous  space  of  characteristics  ji,e.  Replacing  u  with 
e  and  £  with  &  according  to  a  duality  principle  yields  problems 
similar  to  those  discussed  in  VI-A  to  vi-D,  which  can  be  solved  using 
the  functions  developed  in  VI-C  and  VI-D. 
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VII.  SCATTERING  FROM  BODIES  OF  REVOLUTION 


A.  Introduction.  The  objective  of  this  work  is  to  determine  matrix 
methods  of  solving  for  the  surface  currents  induced  on  conducting  bodies 
of  revolution  by  incident  electromagnetic  waves.  These  currents  can  then 
be  us°d  to  calculate  the  scattered  fields. 

X  2 

The  general  problem  has  been  solved  numerically  by  Andreas enXj  for 
bodies  with  maximum  cross-section  perimeters  of  twenty  wavelengths-  His 
met’  >d  consists  of  obtaining  approximate  solutions  of  the  integral 
equations  for  the  induced  currents. 

In  this  section,  the  general  solution  for  the  body  currents  will 
be  formulated  by  the  method  of  moments.  From  this,  two  specific  solutions, 
namely,  a  field  matching  solution  and  a  solution  by  Galerkin *  s  method, 
will  be  described. 


B-  Description  of  Problem.  Consider  an  electromagnetic  wave  incident 
on  the  surface  of  an  arbitrary  conducting  body  of  revolution.  The 
coordinates  for  the  body  and  the  inciden  wave  are  illustrated  in  Fig. 

7-1  and  7-2  respectively.  The  incident  field  (JT )  induces  surface 
currents  on  the  body  which  are  the  sources  of  the  scattered  field  (£  ). 

The  boundary  conditions  require  that  the  total  tangential  £  be  zero  on 


the  body  surface  or 

■(t)  -  *  SjangW 


T?S 

■^tang 


(7-1) 


For  convenience,  (t^,  Zg)  are  defined  as  unit  vectors  of  the  local 


orthogonal  coordinate  system  defined  on  the  body  surface,  t^  corresponds 
to  the  t  vector  and  tp  is  3  in  Fig.  7-1.  la  terms  of  these  coordinates, 
(7-1)  becomes 
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Fig,  7-2  Coordinates  of  incident  Plane  Wove. 


n  =  1.  2 


(7-2) 


tn  *  Is(z)  =  -  \  •  tix) 


The  scattered  fields  are  expressed  in  terras  of  the  surface  currents 


(7-3) 


where  the  operator  L  is  defined  below.  Equation  (T-2)  becomes 


n  *  L  ^(^'5  =  -tn  *  r(r)  n  =  1,  2 


(7-*0 


The  problem  is  now  reduced  to  solving  (7-4)  for  J  (£*)• 

Tc  define  the  operator  L,  (7-3)  is  expressed  in  terras  of  the  magnetic 
vector  and  scalar  potential  s .  The  vector  potential  for  the  induced  surfac 


current  )  is 


’)  -Tf-aa’ 


(7-5) 


body 


where  in  terms  of  the  spherical  coordinates  of  Fig.  7-1 

r  +  r,c  -  2rr ’ (cos  9  cos  9’  + 

.  I1/2 

sin  9  Sin  9'  cos  (0  -  0’))J 


The  scalar  potential  is  given  by 


*S(e)  =  [vf  ‘  J]  ~TT 

body 


(7-6) 


where  *  is  the  differential  operator  tangent  to  the  body  surface.  In 

v 

berms  of  the  local  body  coordinates,  this  operation  on  ar  arbitrary  vector 


£  is 


^t  *  ^  "  3tT  +  Bt  +  f  |fB0 


where  p  is  ths  cylindrical  coordinate  radius  and  v  is  the  counterclockwise 


angle  between  t  and  the  z  axis. 


Now  the  scattered  field  is  expressed  in  terms  of  the  induced  current 


as 


£SU)  =  L  j6T(x)  -  vt  f(x) 


t  i-w 

\  1 } 


where  ’9i.  is  the  gradient  operator  on  the  body  surface 


—  _  f  c  0  d 

Vt  “  °  dr 


D 


The  definition  of  L  is  implied  by  (7-7) 


C.  Expansion  of  the  Incident  Field  and  Surface  Current.  .An  arbitrary 
plane  wave  can  be  expanded  in  a  discrete  set  of  cylindrical  nodes.  Each 
mode  has  a  different  azimuthal  wave  number  and  each  propagates  with  the 
sane  velocity  along  the  z-axis  of  Fig.  7-2.  The  advantage  of  this 
expansion  is  that  the  current  induced  by  each  mode  can  be  computed 
separately  and  the  total  current  can  then  be  found  by  superimposing  all 
of  the  significant  mode  currents.  This  procedure  is  necessary  for 
handling  3arge  bodies  on  available  computers . 

2 

The  details  of  the  expansion  are  given  by  Andreases-  The  propagation 
vector  £  of  the  incident  wave  is  chosen  to  be  in  the  x-z  plane  of  Fig.  7-2. 
Then  the  incident  wave  is  separated  into  a  wave  polarized  parallel  to 
the  x-z  plane  and  a  wave  polarized  perpendicular  to  this  plane. 


The  expansion  for  the  tangential  electric  field  component  of  the 

-»  2 
parallel  polarized  wave  gT  at  the  surface  of  the  body  is 


CO 

Bti(t)  cos  (s0)  +  %  ^ 


T^Ct)  sin  (z0) 


/  __ 


\  i 


T-ii) 


where 


The  total  current  induced  by  this  incident  field  is  then  of  the  form 


.total,  . 

(x) 


[jtl(t)  cos  z4  *  u^Ct)  sin  spj 


n=0  | 


+  0  f'^i(t)  sis  =0  *  JjfcCt)  cos 
L  y-  | 


For  each  incident  rode  the 


neea  dt  eutner  n 


i  I SFS  *S  O— 


Ci(x)  - 1 


t  (t)  cos  np  +  P  u -,{t)  sin 

vi  v-2- 


i 


for  the  parallel  solarized.  case  and 


CoCx)  *  t  JT0(t}  sin  J  --CjCi)  cos  sa0  (7-13) 


roendicular  polarization .  The  rsthods  of  solution  discussed 


bein'*  are  valid  for  solution  for  either  (7-12)  or  (7-15) - 


xt  xs  tsen 


convenient  tc  express  ( 7-12)  and  (7-13)  is  the  general  form 


=  t  oT(t)  h^{30)  +  3  j~(t)  hrfC^p)  (7-1%) 


le  advantage  os  tne  anove  expansions  xs  nov  evident  sxsee  a  sxngxe  setnc-n 


of  solution  can  be  used  to  sens' cutely  solve  for  the  cn 


need  5n 


B.  Solutions  by  Method  of  ?foaents.  She  general  fomailation  for  the 
current  solutions  will  be  obtained  by  the  application  of  the  net-bod  of 
nosen ts  (Section  S-D)  to  (7-1)-  Suitable  choices  of  aggrc-xl  rating  functi 
for  the  cur-iiDt  and  oeiohting  functions  give  the  field  matching  solution 


in,s  vMcb  are  discussed  ^slo». 


The  coordinates  of  the  body  currents  are  called  the  "source  coordinates" 
and  are  the  primed  coordinates  in  Fig*  7-1*  The  coordinates  of  the  incident 
and  scattered  fields  are  denoted  by  the  unprimec'  "field  coordinates." 

From  (7-14),  for  a  given  mode,  say  the  m-th  mode,  the  curre-ci  induced  by 
either  incident  field  polarization  is  of  the  form 

^(r')  =  V  jJ.(f)  ht(m0')  +  0*  f)  h^(m0*)  (7-15) 

A  procedure  analogous  to  the  development  of  (7-1** )  is  used  to  express 
the  incident  field  corresponding  to  (7-15)  &s 

=  t  l^(t)  htfm0)  +  0  E*j(t)  h^(m0)  (7-16) 

where  h_..(in0)  and  h^(m0)  are  the  respective  trigonometric  functions  in 
(7-12)  or  (7-13)  for  whichever  mode  and  polarization  is  being  considered. 

The  scattered  fields  due  to  J 


til)  -  1  ■£(*’)  =  -  ff  4a'  I  k  5-5-  +  V—  Vm 


bo£y 


-jkR  e-jkR 

—  T 


#*■> 


(7-17) 

The  t  and  0  scattered  field  components  are  found  by  taking  the  inner  product 
of  (7-17)  with  the  t  and  0  unit  vectors  at  the  field  points.  The  result 
for  either  polarization  is 


S 


sin  v  sin  v'  (G  .  +  G  -, ) 
vm-l  m+1' 


2  cos  v  cos  v*  GmJ  J^,(t!)  -  sin  v  (G^  -  G^) 

d f+a/ 


a 


m 


•$.  (*’> 


dG_ 


m 


3F"  p 


dt’ 


(7-18) 
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where 


vwid 


?  '  L  #*')  =*  §  [^in  V  <0^  -  Gwl)  J^f) 


p  dt* 


+  (0m-i +  <W  #*'> 

2ra  I"  .  sin  v' 

P 


r  r&jftt')  .  , 


G  p'  dt’> 
nr 


r  e‘ikE° 

G  =  /  -r~=r -  COS  EB0'  * 

m  /  £R  ^  r 

0 


(7-19) 


(7-20) 


Rc  =  R(0  =  0) 


To  apply  the  method  of  moments  to  either  polarization  of  the  m-th 
mode,  the  mode  current  terms  are  approximated  as 

N 

J^(f)  «  Y,  ai  Ti^f) 

i=l 

N 

*  y.  h  $i(t’ }  (T-2i> 

i=l 

vhere  T^(t’)  and  <f>,(t*)are  the  approximating  functions  in  the  method  of 
moments.  A  suitable  set  of  weighting  functions  (w^,  u)^)  and 

an  inner  product  are  defined.  The  approximate  currents  (7-21)  are 
substituted  in  (7-18)  and  (7-19)  and  the  inner  products  of  (7-4)  with 
the  weighting  functions  are  formed  as 
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/J-  *  T  T  /  t  \  I.IV 

r  -  i  i 

si.  r> ■*■/  — .  t  S  I 

j 

4  *  L  jls(£!).  w> 

4  •  w>j 

where  the  current  and  field  terms  are  defined  in  (7-15)  and  (7-1&) -  Both 
of  the  matrices  of  (7-22)  are  complex.  The  left  sides  of  the  respective 
real  and  imaginary  matrices  of  (7-22)  can  be  rearranged  as 


[i^’i 

L  id  J 

‘°i 

♦ 

% 

<t  •  t(x.)>  v 

<t  •  jt U),  y 

h 

4  •  ^U),  ux> 

ti»'] 

*■  -in 

J-ti 

• 

* 

PN 

4  •  /Ce).  V 

i  =  1,2,..., N 

5  -  1,2,.. .,N 

(7-23) 


The  solutions  for  the  complex  currents,  QL  and  ^  are  found  by  inverting 
the  [i]  matrices.  Thus 


(7-24) 


where  the  definitions  of  the  matrices  in  (7-24)  are  implied  from  (7-23)* 


E.  Field  Matching  Solution.  The  field  matching  solution  amounts 
to  satisfying  (7-1)  at  N  points  along  the  body  perimeter  from  t  to 
(see  Fig.  7-1)  for  each  polarization  of  each  mode.  The  perimeter  length 
from  t  to  tb  is  divided  into  N  segments.  The  length  of  the  k-th  segment 
is  At^  and  the  center  point  is  t  as  illustrated  in  Fig.  7-1-  Since  in 
(7-22)  the  0  variations  in  the  scattered  and  incident  fields  are  the  same, 
the  h(m0)  terms  can  be  cancelled  by  choosing  the  weighting  functions  to  be 
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♦  *  *  > 


\  =  6(t  -  t>K)  k  =  1,  2, 


N 


Then  the  inner  product  is  defined  as 


<A,B> 


A  B  dt 


•n 

\  I 


-25) 


t 

a 

where  A  and  B  are  matrices  and  A  is  the  transpose  of  A.  Now  (7-22)  is 
seen  to  be  equivalent  to  (7-1)  at  the  N  points  along  the  body  perimeter. 
The  current  approximating  functions  in  (7-21)  can  be  chosen  to  be 


Tk(f)  =  $k(t’)=  1  on£tk 

0  elsewhere 

This  amounts  to  assuming  that  the  currents  are  of  constant  amplitude  on 
each  interval.  Other  current  expansions,  such  as  a  piecewise  linear 
approximation,  are  possible  but  the  computations  are  more  cumbersome. 

Substitution  of  the  approximate  currents  into  (7-18)  and  (7-19) 
give  approximations  to  the  scattered  fields  at  the  sampling  points. 

These  are  substituted  in  (7-22)  and  the  approximate  current  amplitudes 
are  computed  from  (7-2*). 

The  accuracy  of  the  results  improves  as  the  number  of  sampling 
points  is  increased.  The  scattered  field  at  each  field  point  on  the  body 
is  calculated  as  the  sum  of  the  contributions  to  the  currents  on  each 
ox  the  source  intervals.  For  each  contribution,  (7-l3)  and  (7-19)  must 
be  evaluated.  These  calculations  will  usually  be  approximate.  When 
the  field  point  (t^)  is  several  segments  away  from  the  source  sampling 
point  (t'. ),  sufficient  accuracy  may  be  obtained  if  the  integrands  of 
(7»l8)  and  (7-19)  are  assumed  constant  over  the  sampled  source  intervals 
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and  equal  to  their  respective  values  at  the  center  of  each  source  interval. 
In  general  form,  the  integrals  are  then  approximated  as 

/  I(t')  av  a  Ity 

When  the  field  and  source  segments  are  within  a  few  segments  of  each 
other,  it  may  he  necessary  to  subdivide  each  source  interval  into  sub¬ 
intervals.  The  contributions  to  the  scattered  field  from  the  currents 
on  each  subinterval  can  be  determined  by  evaluating  the  integrands  at 
the  center  of  each  subinterval  and  evaluating  the  respective  integrals 

as  M 

J  **'>  dt'  -  t  I(tw> 

•>=1 

where  M  is  the  number  of  subintervals  of  the  k-tb  source  interval.  When 
t^  =  the  Gfc  terms  (7-20)  have  singularities  since  Rq  =  0.  The 
contribution  from  the  currents  outside  of  the  small  region 

f  -  €  <_  <  0’  +  € 

can  be  determined  as  in  the  above  case.  The  contributions  from  the 
currents  within  the  remaining  small  surface  around  the  field  point  can 
be  determined  by  sin  approximate  analytic  procedure.  This  is  discussed 
in  section  VII-G. 

The  current  derivative  terms  in  (7-18)  and  (7-19)  can  be  approximated 
by  finite  difference  approximations.  For  example,  the  derivative  of  the 
t-component  of  current  at  the  k-th  sampling  point  is 
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6t ' 


Vi  •  °k-; 


2 


Ati 
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Higher  order  approximations  could  be  used  if  necessary. 

The  evaluations  of  fche  G^  terms  defined  in  (7-20)  and  oG^/oi  are 
discussed  in  Section  VII-G. 


F.  Solution  by  Galerkin's  Method.  The  proposed  solution  by  Galerkin's 
method  is  equivalent  to  requiring 

<L  £(£)>  =  -  <gX(r),  ^(x)-  (7-26) 

where  is  the  approximation  to  the  current  on  the  k-th  body  segment. 

To  obtain  this  formulation  from  the  method  of  moments,  the  weighting 

functions  w,  are  equal  to  the  current  approximating  fun-:. cions.  For  example 

let 

\  =  =  \(tf)  =  1  on  At^ 

0  elsewhere. 

The  inner  product  is  defined  by  (7-25).  Equations  (7-22)  are  now  the 
equations  of  Galerkin's  metnod. 

The  constraint  imposed  on  the  solution  by  (7-26)  is  more  stringent 
than  that  of  the  field  matching  solution.  Therefore  it  may  be  possible 
to  obtain  a  higher  degree  of  accuracy  in  the  variational  solution  with 
the  same  number  of  body  intervals  as  in  the  field  matching  method.  However, 
the  calculations  of  the  matrix  terms  in  (7-23)  are  more  cumbersome. 

Many  numerical  techniques  which  are  analogous  to  those  described  in  the 
previous  section  can  be  used  to  simplify  the  computations. 


(j.  EspSSSiC-n  O?  0  **  /p  .  Tn  sithsr  ths  figM  sstcbisg  qt  ysristioDsl 

solution,  one  of  the  major  computations  is  the  evaluation  of  G^  and  bG^/bt 

in  (7~l8)  and  (7-19).  G  is  defined  in  (7-20)  as 

m 


m 


-I 


*  “OkH 


"Sr 


cos 


f*  dy>* 


(7-20) 


0 


where 


r 2  2 

+  r*  -  2rr'  (cos  9  cos  9’  +  sin  0  sin  9*  cos  $*) 


(7-27) 


Using  the  derivative  chain  rule 


BG  cG  \  oG 

m  _  m  or  m  d9 

ou  <3r  Bt  B9-  St 


(7-28) 


Applying  (7-20),  (7-28)  is  expressed  as 


(f,  +  f2  cos  $’)  cos  m0’  d0f. 


ir  -jkRQ 

+  /  — - y- (f^  +  f.  cos  0’)  cos  m0’  d0*  (7-29) 

£  R  ' 

0  o 

where  the  f^*s  are  the  appropriate  functions  of  (r,  r’,  9,  ©’)  which 
result  from  (7-28). 

For  the  integration  techniques  discussed  in  Section  VII-E,  (7-20) 
and  (7-29)  are  integrals  over  the  source  loop  at  the  middle  of  the  source 
interval,  (see  Fig.  7-1)- 

The  integrals  can  now  he  expressed  in  a  general  form  as 


n 


cos  p0'  cos  m0’  d0’ 


m  =  0,  1,  2,  . . . 

P  *  0,  1 
n  =  1,  2,  3 

(7-30) 


when  the  integrands  of  (7-30)  have  singularities  at  Rq  =  0,  approximate 
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analytic  integration  over  the  singularities  is  necessary.  The  details 

p 

of  this  method  are  given  in  the  references.- 

For  the  integrals  which  do  not  include  singularities,  the  complexity 
of  the  integrand  prohibits  exact  analytic  evaluation.  Numerical  integration 
is  possible  but  the  computation  is  cumbersome.  Another  approach  is  to 
expand  e  ^  ^/R"*  in  functions  of  the  integration  variable  0  * .  Substitution 
of  these  expansions  then  allows  approximate  evaluation  of  ( 7-30 )  by 
simple  analytic  integrations.  This  is  the  method  which  is  proposed  for 
the  body  of  revolution  problem. 

To  illustrate  the  expansion,  consider 


R 


n 


cos  kS  .  sin  kR 

-  1 


n 


S 


n 


(7-31) 


Denoting  either  cos  kR  or  sin  kR  as  T(kS),  each  term  in  (7-31)  can  be 
expanded  as 


N 


»  [A  +  3e"°^  +  Ce'#']  V  a^i(0* )  0  <  0’  <  it 

1=0  (7-32) 

for  n  =  1,  2,  3  and  R  >  e  which  is  a  small  positive  constant  to  be 
determined  below.  The  bracketed  terms  of  (7-32)  are  chosen  so  that 


[a  +  Be~°^’  +  Ce“^l  =  -- 
L  i 


0  <  0*  <  « 


(7-33) 


R 


Figure  7-3  is  a  typical  curve  of  1/Rn  versus  0'  for  closely  spaced  loops. 
A,  B  and  a  au*e  chosen  so  that 


-  -O0’  1 

+  Be  ^  -  — 


„n 


>  — 

RU 


at  0’  =  01 >  04,  01 
at  0*  =  0’,  0p 


(7-34) 
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Fig.  7-3  GropWcol  Illustration  of  the  Approximation  of  ^r 
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n 


8(04) 
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B  = 


5(0/,  5 
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B(0p 


n 


~°0U  ~°03 


1  ~°05 

A  =  - -  -  Be  ^ 


b(0:) 


n 


The  constants  €  and  3  are  chosen  so  that 


(7-35) 


^  -  (A  +  Be-0^' )  =  Ce“^ 


and 

Ce  $  — ^ 0  for  0*  >  04 

If  0^  is  properly  chosen,  (?-33)  is  satisfied  when 


C  = 


R(0{) 


n 


-Q0* 

(A  +  Be  X) 


P  - 


— L—  ~  (A  +  Be*^2) 

5(0p)n 


To  complete  the  expansion,  (7-32)  is  rearranged  as 
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(7-36) 


H 


I 


7^1 


A  +  Be 


+  Ce 


«  T(kR) 

Since  the  T(kR)  are  sinusoidal  functions  which  are  even  about  0*  =  s, 
it  is  convenient  to  choose  the  (0  * )  as  trigonometric  functions  with 
the  same  symmetry.  For  example 


constant,  cos  i  0*,  sin  (i  *  i)0’ 

Equation  {  {-3&}  is  a  square  error  approximation  and  the  coefficients 

a.  are  determined  by  satisfying  (7-36)  for  H  values  of  0*. 

Figure  gives  the  results  of  representative  expansions  of  cos  ks/R 
and  cos  kH/R'*  as 

f^(0»)  =  [a  +  Be~°^  -r  Ce~^  j  £aQ  +  a,  cos  05  +  cos  20'  + 


sxn  o_  +  a^  sin^  |1 


(7-375 


Denoting  the  expanded  function  as  f  (0 * ) ,  the  error  at  any  0!^  is  defined 


as 


=(0j!  = 


fM>  *  VW 


max 


From  the  error  curves  in  Fig.  7-4,  the  average  magnitude  of  error  for 
either  expansion  is  less  than  0-3  percent.  For  these  cases,  E  = 

*  lulU 

O.OoTX.  It  was  found  that  for  R_.„  >  O.067 X,  the  average  error  decreased. 


mm 


so  the  cases  presented  give  an  estimate  of  the  maximum  errors  that  can 
be  exoected  for  expansions  of  bodies  of  revolution  with  diameters  less 


Fig.  7-4  Results  of  Computed  Expansions 


than  O.p^  and-  for  B  .  >  €  =  0.067*-.  For  larger  diameters,  it  may  be 

mm  — 

necessary  to  increase  the  number  of  terms  in  (7-37) • 
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VIII.  inversion  of  matrices 


A.  Introduction.  The  inversion  of  a  matrix  is  closely  related  to 
the  solution  of  a  system  of  linear  algebraic  equations.  The  n  lineal 
equations  with  n  unknowns 


allxl  +  a12x2  +  *  *  ‘  +  alnxn  ~  yl 


a21Xl  +  a22X2  +  ' 


‘  +  a2nXn  =  y2 


(8-1) 


a  ,x.  +  a  ,_x0  +  .  .  .  +  a  x  =  y 
nl  1  n2  2  nn  n  Jn 


may  be  written  as  a  matrix  equation 

Ax  =  y  (8-2) 

where  A  is  an  n  x  n  matrix,  x  =  (x^,  xp,  ...,  x„)  an  unknown  vector 

and  y=  (y  ,  y  ,  . ..,  y  )  a  known  vector.  This  matrix  equation  can 
JL  c.  21 

be  conveniently  written  as  a  partitioned  matrix  (Ajl)  with  a  column 
vector  [-— ] 

**y 

(A|t)  [2L]  =  0  (8-5) 

V 

where  I  denotes  an  identity  matrix. 

The  elf.  nentary  row  operations  are  defined  as 

a)  multiply  any  row  by  a  constant) 

b)  interchange  any  two  rows; 

c)  add  one  row  multiplied  by  a  constant  c  to  any  other  row, 
or  their  combination. 

When  the  elementary  row  operations  are  performed  on  (8-1),  or  equivalently, 
on  the  partitioned  matrix  ( A j I)  of  (8-3)  they  do  not  change  the  solution 
of  the  system  of  Eqs.  (3-1)  or  the  matrix  equation  (8-3)* 
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Let  L  be  a  matrix  representing  some  combination  or  product  of  the 
elementary  operations  described  above.  After  these  operations  on  (8-3), 
we  have 


(la|l)  [J]  =  o 

V 


(8-4) 


If  we  carry  out  these  elementary  row  operations  so  that  LA  becomes  an 
identity  matrix,  i.e.,  LA  =  I,  then  the  inverted  matrix  A  1  is  found  to 
be 


A"1  =  L 


(8-5) 


B.  Method  of  Gauss-Jordan  Reduction.  The  method  of  Gauss-Jordan 
reduction  is  a  systematic  procedure  designed  to  execute  the  method  outlined 
in  Section  A.  Let  A  be  an  n  x  n  matrix  to  be  inverted.  The  partitioned 
matrix  of  (8-3)  can  be  written  as  an  n  x  2n  matrix 


*11 

a12 

«  •  •  -» 
l,n-l 

aln 

1 

0 

0 

0 

a21 

OJ 

af'J 

*  *  *  a2,n-i 

a2n 

0 

1 

0 

0 

• 
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• 

• 

• 

• 

• 

• 

• 

* 

* 

♦ 

• 

• 

• 

4 

♦ 

* 

• 

• 

• 

• 

» 

• 

Vi,i 

an-l,2 

•  *  *  Bi  m.  , 

n-x,n-l 

a  T 

n-l,n 

0 

0 

1 

0 

anl 

an2 

•  *  4  2*  - 

n,n  -x 

a 

n,n 

0 

0 

•  44  0 

I 

(8-6) 


Search  the  ^irst  column  a.,  (i  =  1,  2,  n)  for  the  largest 

element  in  magnitude,  say 

\\l\  -  laiil  ^  ~  2>  **•»  n)  (8-7) 

If  =  0,  then  the  first  column  of  A  is  zero,  A  is  a  singular  matrix 
and  hence  A  has  no  inverse  in  the  usual  sense.  Therefore  we  may  assume 
a^-,  ■£  0.  Interchange  the  first  and  k-th  rows,  and  then  divide  the  firs' 
row  through  by  a^.  If  the  interchange  of  the  first  and  k-th  rows  is 
followed  by  an  interchange  of  the  first  and  k-th  columns  of  the  second 
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submatrix,  then  the  resulting  matrix  becomes  more  orderly.  To  make 


sure  that  the  interchange  of  columns  of  the  second  submatrix  does  not 
disturb  the  matrix  equation  (8-3)  or  (8-k),  it  is  required  that  and 

be  interchanged  also.  Finally  we  use  the  first  row  to  eliminate  all  but 
the  first  term  of  the  first  column  by  the  elementary  row  operations.  Thus 


(L(1)A|L(1))  [-JL-]  ,  o 


(8-8) 


where  'y  [M^  is  a  permittation  matrix 

which  transforms  y  =  (y1 ,  y?,  y  )  into  y^  =  (y,„.  yp}  y,,  ... 


'!>  2 
\  i 


yk-l’  vl>  yk+l>  Vj 


1 

b12 

°l,n-l 

bln 

0 

...  0 

0 

0 

• 

* 

°22 

• 

• 

b2,n-l 

♦ 

* 

b2n 

• 

* 

b21 

• 

* 

1 

• 

• 

...  c 

• 

• 

0 

• 

• 

• 

0 

• 

d 

n-1,2 

• 

Dn-l,n-l 

• 

b  i 
n-l,n 

• 

Dn-l,l 

• 

0 

♦ 

...  1 

♦ 

0 

0 

b  _  ... 

n2 

b 

n,n-l 

b 

n,n 

bn,3 

0 

*  * «  0 

1 

(8-9) 

The  elements  of  the  matrix  (8-9)  can  be  written  in  terms  of  the  elements 
of  the  matiix  (8-6)  by  inspection. 

We  proceed  to  obtain  a  new  mat:  ix  equation 

(L(2)a|L(2))  C-jgy]  =  0  (3-10) 

(o\ 

from  (8-8)  by  the  elementary  row  operations  such  that  the  submatrix  Lv  'A 
has  zero  elements  for  the  first  and  second  column  except  the  first  and 
second  diagonal  terms  which  are  normalized  to  unity.  This  can  be  achieved 
by  the  following  steps j  : 

1)  Find  an  integer  k  such  that 

1*^2 1  ^  l^ig !  =  2;  5}  •••}  n-1.  n) 
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The  lb,  0|  may  be  assumed  to  be  non-zero  otherwise  the  original  matrix 


is  singular. 

2)  interchange  the  second  and  k-th  rows  and  also  interchange  the 
second  and  k-th  columns  of  the  second  submatrix * 

3)  divide  the  (new)  second  row  by  b^, 

4)  use  the  second  row  to  eliminate  the  second  column  by  the 
elementary  operations  except,  of  course,  the  second  element  of  the  column, 

5)  interchange  y^  and  y  . 

It  is  readily  seen  that  these  steps  lead  to  the  matrix  equation  (8-IO). 
Assuming  that  the  matrix  A  is  nonsingular,  then  after  repeating  the 


process  n  times,  we  obtain 


(L(n)AiL(n))  [-— r]  =  0 


(3-11) 


where  l/n^A  is  an  identity  matrix,  and  is  a  vector  resulting  from 
the  rearrangement  of  the  original  y.  Let  M  be  a  permutation  matrix 
transforming  y  into 


=  My 


(3 -12) 


(n)  ~(n) 

Since  M  is  a  permutation  matrix  so  that  MM  =  I,  L  and  L  are  related 


L(n)  _  £(n)?4 


(8-13) 


Equations  (8-12)  and  (8-13)  imply  that  if  y.  and  y..  are  interchanged  in 
(n)  -  ^rn\ 

,  then  b  and  b  of  Lv  1  (k  =  1,  2,  n)  must  be  interchanged. 

Kl  KJ 


C.  Commutative  Matrices.  When  a  boundary  value  problem  has  some 
degree  of  symmetry  such  as  rotational  symmetry  or  reflectiona  1  symmetry, 
the  matrix  obtained  by  the  method  of  approximation  described  in  previous 
chapters  often  exhibits  a  special  structure.  This  special  structure  of 
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the  mati  ix  can  usually  be  characterized  by  the  fact  that  this  matrix 
commutes  with  another  matrix  which  represents  the  nature  of  the  symmetry. 

The  product  of  two  matrices  A  and  B  depends ,  in  general,  on  the  order 
of  the  factors,  i.e.,  AB  ^  BA.  In  the  very  special  situation  when  the 
product  is  independent  of  the  order  of  the  factor  i.e.,  AB  =  BA,  t  en 
the  matrices  A  and  B  are  said  to  commute.  A  case  of  special  interest 
to  us  is  when  a  given  matrix  commutes  with  a  diagonal  matrix. 

Any  matrix  that  commutes  with  a.  diagonal  matrix  having  distinct 
diagonal  terms  is  necessarily  a  diagonal  matrix. 

To  prove  this  statement  let  A  =  [*2.  -1  be  an  arbitrary  n  x  n  matrix 

3-d 

and  D  be  an  n  x  n  diagonal  matrix  with  distinct  diagonal  terms; 


D  =  [d.  6.  .],  d.  £  d. 
i  ij  i  r  Z 


it  it  6 


Since  AD  =  DA,  we  have 


[a.  J[d.6. .]  =  [a.8.  .][a,  .3 

L  13  i  xo  x  io  J.0 


which  is  equivalent  to 


“ik  8kj  V  ■ 


n 

di  6ik 

k=l 


This  equation  leads  to 


[a.^(d  -  d  )]  =  [0] 

— vi  u  — 

where  [o]  denotes  a  zero  matrix.  The  above  equality  holds  only  if 


a.  .  =  0 
10 


if  d.  t  d.  for  i  t  3 

i  o 


which  shows  that  the  matrix  A  =  [cs..3  must  be  diagonal. 

It  follows  from  the  proof  given  above  that  if  a  matrix  A  commutes, 


with  a  diagonal  matrix  D  whose  diagonal  terms  of  equal  value  appear 
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o/t.%'%4*  -?  ■*  ro  1 

W :»uCuu.ur  t  VJ.J  y 

then 

Tnof  t«^  v 

mu.£  v  is  w 

2ri.x  wi_tn  zsi*o  off  * 

✓* 

4  r .  1  n 

diagonal  suDmatrice 

s-  ''or  example,  ii' 

A  is  a  0 

x  0  matrix  ia.  .  j  ana  u 

is  a  diagonal 

matrix; 

! 

r 

-1 

D  = 

dl 

0 

0  0 

0 

0 

0 

d2 

0  0 

0 

0 

0 

0 

d  0 

0 

0 

0 

0 

u 

0  d5 

0 

0 

(8- 15) 
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0 

0  0 

d, 

y 

0 

0 

0 

0  0 

0 

d5 

then  AD  ■=  DA  implies  that 

A  = 

“ll 

0 

0 

0 

( 

3 

0 

0 

a22 

a23 

0 

{ 

3 

0 

0 

°32 

c 
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3 

0 
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°46 
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^6 

yy 

0 

— 

0 

0 

a61 

o5 

a66_ 

It  is  well  known  that  a  matrix  equation  remains  invariant  under  a 
similarity  transformation  which  includes  unitary  and  orthogonal  trans¬ 
formations  as  its  special  cases.  For  example  the  matrix  equation 


AC  =  CA  (5-15) 

which  denotes  the  commutahiiity  of  the  matrices  A  and  C.  after  a 
similarity  transformation  (by  a  non-singular  matrix  P)  takes  the  same 
form 

EB  =  DB  (S-l6) 

where 
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and 


(3-17) 


B  =  P~J'AP 


D  =  P_1CP 


Now  suppose  that  C  of  (8-15)  can  be  diagonalized  by  a  similarity 
transformation  and  that  the  diagonalized  matrix  D  has  distinct  diagonal 
terms-  Then  the  matrix  B  of  (8-l6)  obtained  by  the  same  similarity 
transformation  from  matrix  A  must  be  diagonal  also.  On  the  other  hand, 
if  D  is  a  diagonal  matrix  with  some  repeated  diagonal  terms  as  shown  in 
(8-13),  then  B  must  be  a  supermatrix  with  zero  off-diagonal  subsatrices 
as  shown  in  (8-14). 

The  invariant  property  of  a  matrix  equation  under  the  similarity 
transformation  sometimes  can  be  used  to  simplify  the  inversion  of  a  matrix. 
Let  A  be  a  matrix  to  be  inverted.  Suppose  we  can  find  a  matrix  C  which 
commutes  with  A  and  whose  eigenvalues  and  the  corresponding  eigenelements 
can  oe  easily  found  so  that  a  nonsingular  matrix  P,  which  will  diagonalize 
C  by  a  similarity  transformation,  can  be  obtained.  In  addition,  if  the 
diagonalized  matrix  has  distinct  diagonal  terms,  then  the  matrix  A  is 
reduced  by  the  same  similarity  transformation  to  a  diagonal  matrix  which 
can  be  inverted  easily.  The  desired  inverted  matrix  A  1  can  be  obtained 
by  reversing  the  original  similarity  transformation.  In  case  the  diagonal 
terms  of  the  diagonalized  matrix  obtained  by  the  similarity  transformation 
are  not  distinct,  the  matrix  A  is  reduced  by  the  same  similarity  trans¬ 
formation  to  a  supermatrix  with  zero  off-diagonal  submatrices.  The  inverted 
matrix  A  *"  is  then  obtained  by  first  inverting  ~hese  diagonal  submatrices 
and  then  reversing  the  original  similarity  transformation. 

The  practicability  of  the  method  described  above  depends  on  whether 
we  can  find  a  second  matrix  which  commutes  with  the  matrix  to  be  inverted 
and  whose  eigenvalues  and  eigenelements  are  known  or  easily  obtainable 
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so  that  a  nonsingular  matrix  can  be  constructed  to  diagonalize  the  matrix. 


To  find  such  matrix,  we  consider  the  specific  configuration  of  a  four- 
pert  network  shown  in  Fig.  8-1.  It  is  easily  seen  that  when  the  network 
undergoes  a  rotation  of  s/2,  s,  or  3*/2,  the  electrical  structure  remains 
unchanged.  Let  e  =  (e1,  e^,  e,,  eSi)  denote  the  set  of  excitations,  r 
the  set  of  corresponding  responses  in  the  network  due  to  the  excitation 
and  S  the  matrix  representation  of  the  network  parameters  determined  by 
the  electrical  structure  of  the  network.  Then  the  equation  relating 
the  excitations  and  responses  of  the  network  is  given  by 

r  =  Se  (8-18) 

Let  R  be  a  rotating  operator  signifying  a  rotation  of  s/2  in  the 
counterclockwise  direction-  A  rotation  of  the  s^t  of  excitations  in  the 
counterclockwise  direction  by  an  angle  of  s/2,  keeping  the  network 
stationary,  gives  a  set  of  new  responses 


r’  =  SRe 


(8-19) 


On  the  other  hand,  the  set  of  new  responses  r*  is  equal  to  the  set  oi 
responses  obtained  by  rotating  the  whole  system,  including  both  the 
excitation  and  network,  counterclockwise  by  an  angle  of  s/2; 


r '  =  Rr  =  RSe 


Thus  by  (8-19)  aud  (8-20),  we  find,  for  ev< 


er/  excitation  e 


Swo  = 


wmen  Ag  fv-sit, 


SR  =  RS 


C  8-21'* 
\ w  t-i-/ 


LO-«  } 


Usually,  the  S  matrix  which  represents  the  electrical  structure  of 
the  network  or  the  dynamic  properties  of  a  boundary  value  problem  is 


rather  complicated  while  the  E  matrix  which  characterises  the  symmetry 
of  the  network  or  the  problem  is  relatively  simple.  This  feature  is 
especially  evident  when  the  problem  has  a  high  degree  of  symmetry. 


D.  Example.  Circular  Loop  Antenna.  The  boundary  value  problem  of 
a  circular  loop  antenna  possesses  a  high  degree  of  symmetry,  namely 
the  rotational  symmetry  about  the  center  of  the  loop.  One  of  the  so^-t 
effective  approaches  for  obtaining  an  approximate  solution  is  to  divide 
the  circular  loop  in  N  equal  sections  and  regard  the  antenna  as  an  n-port 
network.  The  mafcrzx  equation  obtained  by  approximating  the  Integra- 
differential  equation  describing  the  boundary  value  problem  is  a  n  x  n 
matrix  (impedance  matrix) 

A  =  [a  ]  (8-25) 

+3 


The  problem  is  to  invert  the  matrix  A  to  obtain  an  admittance  matrix. 

It  is  seen  from  Fig.  3-2  that  the  boundary  value  problem  is 
invariant  with  respect  to  a  rotation  of  the  system  by  an  angle  2sk/ll 
(k  =  1,  2,  — ,  II)  in  counterclockwise  sense.  Since  this  rotation  can 
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Next  we  consider  the  diagonal! cation  of  the  C  matrix .  For  this 
jurpose  it  is  necessary  to  obtain  the  eigenvalues  and  eigenelements  of 


C  matrix.  The  characteristic  equation 


(C  -  Xl)x  =  0 


(8-2 6) 


gives  the  N  eigenvalues  and  the  corresponding  eigenelements 


.  2nlr 
^  N 

e  ,  k  =  0,  1,  2,  N- 1 


(8-27) 
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L-  zr 


.  2i±£ 


eA>  k-0,  1,  2,  N-l | 


(8-28) 


whv re  e^  denotes  a  unit  element  in  the  N-dimensional  Euclidean  space. 

The  diagonalizing  matrix  P  can  be  obtained  from  the  set  of  normalized 
eigenelements  by  regarding  them  as  the  column  elements  of  P: 


P  =  [7W]  =  -^r  [e  N  3  (M  =  0,  1,  2,  N-l) 


(8-29) 


The  inverse  of  P  is  seen  to  be 

-•  2rtk£ 

P_1  =  P*  =  J_  [e  °  N  ; 


P  =  P*  =  -  [e  W  ]  (k,£  =0,  1,  2,  N-l) 

W 

(8-30) 

Tl.^t  Pisa  unitary  matrix  can  be  seen  from  the  following  relation; 
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The  mat?  x  C  can  be  diagonalized  by  the  following  similarity  transformation 


j  — 

P-1CP  =  [e  N  &ki3NxN  (k,i  =  0,  1,  2,  N-l) 


(8-3D 
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which  has  N  distinct  diagonal  terms  equal  to  the  eigenvalues  of  C.  It 
follows  that  B  =  P  1AP  must  be  a  diagonal  matrix  also. 

We  now  evaluate  the  diagonal  terms  of  B.  The  elements  of  B,  in  general, 
can  be  calculated  by  the  formula 
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(8-32) 


According  to  the  results  in  Section  C,  we  see  that  the  off  diagonal 
terms  vanish,  i.e., 

=  o  ifi/« 


For  the  diagonal  terms,  we  have 


N-l  .  2rt(k-i)l 

ht  -  v  I  “ik  e 

i  ,k=0 

This  expression  can  be  simplified  by  making  use  of  the  commutability  of 
matrices  A  and  C. 

It  follows  from  the  relation  AC  =  CA  that  the  elements  a.,  of  A 

is: 

depend  only  on  the  relative  positions  of  i  and  k.  Let  p  denote  the 
relative  position  of  i  and  k  and  write 


a 


P 


(8-3*0 


A  moment  of  careful  examination  of  (8-33)  and  Fig.  8-2  reveals  that  (8-33) 
may  be  regrouped  as  follows; 


-l6p- 


Integral  number-  assumed 
bv  n 


*  t|] 


The  partial  sum  of  (8-33)  corresponding 


+  n  Y\ 


.  2jt^p 

n  J  ~r 

?aP  e 


.  2rtip 


N  Ct  e 
P 


if  K  is  even 


if  N  is  odd 


[§]  +  i 


.  2it£p 


K  a  e 
P 


.  2*ip 


N  a  e 
P 


-N-.  i 
12!  -  1 


N  a  e 
P 


-  2jt-gP 
*  N 


.  2«^p 


?  ap  6 


if  N  is  even 


Ha  e 
P 


if  N  is  odd 


Where  the  square  bracket  [Q]  denotes  the  maximum  integer  of  Q,  e.g., 

[2  4 13-^  ]  =  2. 


Summation  of  these  terms  over  p  yields 
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(8-35) 


If,  in  addition,  the  reciprocity  property  holds  as  in  the  problem 
of  circular  loop  antenna,  then 


Qik  =  °ki 


(8-36) 
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Consequently ,  (8 -3p)  can  be  written  as 


[N/2] 


2  -  6(p) 


Nvi  2jtin 

— }  >  a  cos  — — =- 
2'j  p  N' 


(8-38) 


which  represent  the  diagonal  terms  of  matrix  B.  Thus  after  the  similarity 
transformation  of  matrix  A,  we  obtain  a  diagonal  matrix 


B  =  P~  A?  =  Sj 


(8-39) 


The  inverse  of  this  matrix  is  easily  seen  to  be 


b"1  =  p'1  a"1  p  =  [,—] 

hi 


Therefore,  A  1  is  given  by 
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Thus  the  elements  of  A  can  be  identified  as 
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(8-62) 


It  suffices  to  evaluate  the  elements  of  the  first  row  of  A  ,  since  the 
elements  of  the  other  rows  can  be  constructed  from  these  elements.  For 


the  first  row,  we  have 

N-l 
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(8-43) 
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Again  for  those  problems  which  satisfy  the  reciprocity  property,  ye  have 
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(8-4 4) 


where  are  given  by  (8-38} . 

As  an  example,  we  consider  the  inversion  of  a  4  x  4  symmetric 


matrix 
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which  commutes  with  a  rotation  matrix  operator 
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By  (8-38) ,  matrix  A  can  be  diagonalized  by  a  similarity  trans¬ 
formation  to  give  a  diagonal  matrix  3  with  the  following  diagonal  terms; 
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44)  the  elements  of  the  first  row  are  found  to  be 
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The  inverted  matrix  A  is  found  according  to  (S-1p)  with  the  elements 
of  the  first  row  given  by  (8-48) ; 
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E.  Example.  Linear  Antenna.  As  another  application,  we  consider 
the  problem  of  a  linear  antenna.  The  electrical  property  of  the  antenna 
remains  unchanged  by  rotating  the  antenna  with  respect  to  the  mid-point 
by  an  angle  rr.  The  matrix  equation  intended  to  approximate  the  integro- 
aifferential  equation  describing  the  boundary  value  problem  can  be  arranged 
by  a  proper  choice  of  the  approximating  functions  (elements)  to  have  a 
matrix  A  =  [<X  .  ]  which  commutes  with  the  matrix  representing  the  rotation; 


that  i s 


AR  =  RA 


(8-50) 


(8-51) 


We  shall  show  that  the  inversion  of  matrix  A  can  be  reduced  to 


the  inversion  of  its  diagonru.  rubmatriees.  To  be  specific  let  A^, 
and  I^n  denote  matrices  of  order  2n.  The  characteristic  equation  for 


*2n 


(8-52) 
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gives  the  eigenvalues  +1  and  -1.  each  of  which  repeats  n- times. 


For  the  eigenvalue  X  =  1  and  x  *  (g^,  f^,  . g^,  . ..,  g^),  tiie 
equation 

(S2n  -  XI2n)x  "  0  <8-53> 

yields  the  following  relations  to  be  satisfied  by  the  eigenelements ; 

El  '  *2n  "  0 

52  '  S2n-1  =  °  (8-54) 

^n  ”  ^n+1  “  ® 

The  corresponding  eigenelements ,  after  normalization,  can  be  chosen  as 
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Similarly  for  the  eigenvalue  X  =  -1,  we  find  the  n  normalized  eigenelements 
as  follows; 
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The  matrix  Ppn  which  diagonalizes  by  a  similarity  transformation 

may  be  cot  structed  from  these  eigenelements  by  regarding  then  as  column 

elements  of  P_  ; 

2n 
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